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ABSTRACT 


The Low Power Atrospheric Compensation Experiment (LACZ) 
spacecraft was launched for NRL in february 1990. The LACE 
flight dynamics experiment will provide on-orbit systems 
identification of the LACE spacecraft. The experiment is 
designed to measure modal frequencies, damping ratios, and 
oscillation amplitudes of the LACE spacecrait. The purpose of 
this study is to develop a finite element model of the LACE 
spacecraft and conduct a dynamics analysis to determine 
natural frequencies and mode shapes. Four configurations of 
the spacecraft are analyzed. This data will be compared with 
actual orbital data and will provide an opportunity for 
improvements in the accuracy of computer simulations of 
flexible structures and multi-body dynamics. Thermoelastic 
effects due to differential heating are addressed to check the 
magnitude of deformations that may cause a problem for 
stability or on-orbit identification. The final phase of this 
study is to conduct a parametric analysis of the spacecraft 
boom to investigate the presence of chaotic vibration for 


combinations of excitation amplitude and frequency. 
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I. INTRODUCTION 


A. PURPOSE 

This study is concerned with the modeling and analysis of 
the Low Power Atmospheric Compensation Experiment (LACE) 
spacecraft. One of the missions of the LACE spacecraft is to 
conduct and obtain flight data for on-orbit system 
identification. The flight dynamics experiment is designed to 
measure modal frequencies, damping ratios, oscillaticn 
anplitude of the LACE spacecraft and vibration intensity 
generated py boom deployments and retractions. This experiment 
will provide an opportunity for improvements in the accuracy 
of computer simulations of large flexible space structures and 
Malti-body dynamics. 

It will also provide a mechanism for evaluating influence 
of magnetic torques, gravity gradient torques and atmospheric 
drag on the LACE-type structures. 

The purpose of this study is to develop a finite element 
model of the LACE spacecraft using the finite element program 
Graphical Interactive Element Total System (GIFTS). Dynamic 
analysis is performed on the model to determine the natural 
frequencies and mode shapes. Thermoelastic effects due to 


differential heating are addressed to check if the magnitude 


of deformations could cause a problem for stability or system 





identification. Finally, a parametric analysis is conducted on 
a model of the spacecraft boom to investigate the possibility 


of chaotic vibrations that may be induced during the mission. 


B. OVERVIEW/BACKGROUND 
1. Spacecraft Description 

NRL-developed LACE was successfully launched on February 
14, 1990, from Cape Canaverai on a DELTA II launch vehicle. 
The spacecraft incorporates three deployable/retractable booms 
of maximum length 150 feet. The 2,800 1b. LACE spacecraft is 
stabilized by a 150 ft. zenith directed gravity gradient boom 
mounted on top of the spacecraft, a momentum wheel with axis 
along the pitch axis and a magnetic damper at the tip of the 
gravity gradient boom. The retroreflector boom is mounted 
forward and deployed along the velocity vector, while the 
balance boom is mounted and pointed aft. Figure 1 shows the 
basic configuration of the spacecraft. 

2. Experimental Hardware 

The flight dynamics experiment hardware consists of 
three germanium corner cubes mounted on the lead boom, on the 
bottom of the bus and on the aft balance boom to serve as 
targets for the 10.6 micron Firepond laser radar at MIT 
Lincoln Laboratory, Lexington, Massachusetts. The Firepond 
laser has a 4 millisecond square wave pulse at a frequency of 


62.5 Hz and pulse energy of 3.2 joules. The Firepond laser 


radar will illuminate the cubes to measure the relative motion 
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Figure 1. LACE Spacecraft Configuration 

of the boom with respect to the main body. Reflections from 
the corner cubes wi’l give differential Doppler information on 
the magnitude of the displacement rates due to boom flexure. 
(Ref. 1) 

The Ultra-Violet Plume Instrument (UVPI) will be used 
to measure the absolute bus rotation rate. It has the 
capability of resolving angular velocities of 5 * 10°5 
radian/sec. (Ref. 2] 


Figure 2 depicts the LACE flight dynamics experiment. 








Figure 2. LACE Flight Dynamics Experiment 


3. Acquisition Sequence 

The recommended acquisition sequence for the LACE 
spacecraft calls for four different boom deployments. 
Initially, the gravity gradient will be extended to 75 feet. 
After the spacecraft has stabilized, the boom will be deployed 
to 150 feet and the momentum wheel will be spun up. Then the 
leading (retroreflector) and tracking (balance) booms will be 
extended to 119.5 feet. The final configuration will have the 
retroreflector boom and gravity gradient boom extended to 150 


feet and the trailing balance boom extended to 75 feet. 





ee me me 


4. Orbital Parameters 
Table = lists the orbital parameters for LsCz. 


Table I. LACE ORBITAL PARAMETERS 


95.47 min 


43.087 decrees 
-00208 


Sezi-major Axis | 6939.351 kn 





C. MOTIVATION 





During the past five years, systen 
flexible space structures has eserged as an important problen. 
Many proposed space missions will involve large space 
structures that are very ‘flexible, contain thousands of 
structural elements and have special mission reauirements, 
such as pointing accuracies. Some of these structures include 
space defense platforms, solar power stations and manned 
laboratories, such as the space station. These structures may 
require some type of active control to carry out ongoing 
maneuvers, suppress and control vibration and achieve accurate 
and reliable pointing. An obstacle to meeting some of these 
objectives can be attributed to the inability to analytically 
model the structural dynamics of these highly flexible 
structures with a high degree of confidence or precision. 


Typically, such information is also required in the design of 


vibration suppression and control systems. 





The inability cto xodel dynamics of iaerce Flexidle 

structures is caused by four Zactors: 

1. Because of leunch costs, lexce space structures are 
constructed ort of very licht cox=posites and cannot 
support their weight in cravity. This prevents 
Ground-based testinc. 


2. New composite materials will vesult in modeling 
uncertainty of the material properties. 


3. Complications arise from the cozbination of very high 
flexibility and large thermal grauients. 


4. On-orbit structure will be exposed to harsh space 
environzents which include particle radiation, solar 
eiiects, gravitational anomalies, extreme 
temperatures and near vacuum. The structure may 
underco physical parameter chances. [Rez. 3] 

General purpose structural modeling and muiti-body dynamics 
computer programs such as GIFTS and TREETOPS, respectively, 
can be used to define and investigate complex space 
structures. However, there is no assurance that optinum 
structural models can be generated with these programs. It is 
envisioned that the on-orbit data received from LACE can be 


used to improve these computer models. 


D. THESIS OUTLINE 

This study contains eight chapters. The second chapter 
provides a discussion of the theoretical basis adopted by 
GIFTS to develop a finite element model and provide dynamic 
analysis. Chapter III describes the development of the simple 


beam model and complex finite element model of the LACE 


spacecraft. Chapter IV provides the natural frequencies and 
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mode shapes cf the spacecraft. in Chapter V a thermoelastic 
analysis is conducted on the LACE spacecraft booa to deternine 
the effects of differential heating. Chapter Vi discusses 
multi-body dynazics and describes a multi-body progran called 
TREETOPS. Chapter VII gives a description and a method to 
identify chaotic vibrations. A parametric analysis of the LACE 
spacecraft boom modeled as a single degree of freedom system 
is provided to determine if chaotic vibrations may be induced. 
The model uses experzmental data for stiffness, mass and 
dGaxping and examines the behuvior for possible chaos in the 


system. Chapter VIII contains conclusions and scope for future 


research. 








II. THEORETICAL FORMULATION 


&. FINITE ELEMENT METHOD 

The finite element method is a numerical procedure to 
compute the response of complex structures. The basic unit of 
this analysis, the discrete finite element, is a geometrically 
sinplified representation of a small part of the physical 
structure. The finite element method views the complete 
complex structure as an assembly of a finite number of these 
finite discrete elements (beams, rods, plates, etc.), each of 
whose properties and deformation responses are Simple, as 
compared to the complete structure. The division of the 
discrete elements is natural, in general, and follows the way 
the actual structure is built (truss members, frames, etc.). 
The elements interconnect at nodes where the elements meet and 
move in unison only after compatibility requirements are met. 
This assures that adjacent elements will not overlap or 
separate. 

Each node has six degrees of freedom with three deflections 
and three rotations. Based on how the elements are connected 
at the rods, the computed properties of the individual 
elements are determined and assembled to obtain the 


equivalent, but more complex properties of the entire 


assembled model. The structural model created can then be used 
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to predict the behavior of the real structure. In this study, 
the structural model will be used to determine the structure’s 
natural Zreauencies and mode shapes. The finite element types 
used to model the LACE spacecraft include rod, beam and 
plate/sheli elements. These elements will be discussed further 
in accordance with the theory adopted used by GIFTS. 
1. Rod Elements 
Rod elements are pin-jointed truss members that are 


capable of carrying axial forces only. Figure 3 shows a 


typical rod element. 





Figure 5. Rod Elements 


The element 1s located in tne (x~y-z) coordinate system, which 
is also known as the reference or global coordinate system. 
Ccordinates along the element are the local or element 
coordinate system (x-Y-2). 

The following derivations can be found in more detail 
in Allen and Haisler [Ref. 4: section 7.2]. The element 
stiffness matrix, derived from the principle of virtual work, 


in local coordinate system, is defined by the equation 
Y y 





paeny 


q - BEI 7 (1) 
i = 28[7 2] 


The matrix [K] represents the truss element stiffness matrix 


in the local coordinate systen. 


Figure 4 shows the rod element oriented at a positive 


angie 8 relative to the global axis. 





Figure 4. Rod Element at 6 


Element and global displacements at nodes 1 and 2 are related 
by 


uy 
4) _feos@ sin@ 0 O |} Ye (2) 
af | o 0 cos@ sind] ju, 

uy, 


or 





{a} = (7) {u} (3) 


ae tm OF 
nr nent 
5 
ete een 
: 


where 


cosé sind 0 0 (4) a ‘ 


{7} = 0 0 cos6 sin® , ‘ 


Using the strain energy equation and the matrix product 


transpose rule, the following relations are derived: 


u= $0 e580) 
= S (ul [K} {u} (5) 
where \ 
[K) = (TITIX (7 (6) 


Therefore, using equations (1) and (6), [K] is transformed to 


the global coordinate system to obtain 


c? cs -c* -cs 
_ Ag| CS S? -cs -s? (7) \ 


Ks . 
Li-c? -cs c? es | ; 

-cs -s? cs 3? | 

where c = cos@ ands = sin6é. 


2. Beam Elements 
Beam elements are characterized as members that are 
capable of resisting bending. The beam element will carry 


axial forces, shear forces and bending moments. Figure 5 


depicts a typical planar beam finite element. } 








Figure S. Beam Element 


The following can be found in more detail in Allen and 


Haisler (Ref. 4: section 7.3.]} 


The stiffness matrix 





EA 

7 0 

0 12ETI 
L* 

0 e 

(Ki =] _ ps 

L 

0 -12ET 
L? 
§2I 

0 =F 


where I = Ty. 


In order to assemble the stiffness matrices, they must 


be transformed to global 


transformation, as before, 





is defined by the equation 





-EA 
a 0 
L 0 
6EI 9 x12EE SEI 
L? L L? 
4EI -6EI 2ET 
L x L 
EA 
0 = 0 0 
L 
-SEI 4 12EI ~-6ETI 
L? LS L? 
2EI 4 -6EI 4AET 
L L? L 


coordinates. The local-global 


is given by 


(8) 


nt ob 


eee 


001000 (9) 


where c = cosd ands = sin®é, 

The remainder of the finite element formulation is 
similar to the previous section on rod elements. The stiffness 
matrix in global coordinates is constructed using equations 
(6), (8) and (9). 

3. Plate Bending Theory 
A flat plate, like a beam, supports transverse loads and 


offers resistance to bending. Figure 6 shows stresses that act 


on a homogeneous linearly elastic plate. 





Figure 6. Plate Stresses 


The normal stresses o, and o, vary linearly with 2 and 


contribute to bending moments M, and M,. 


13 
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The normal stress 0, is negligible when compared with 


o,, 6, andt,,. The transverse shear Stresses t,, andt,. vary 


y 
quadratically with 2. Plate bending in this analysis refers to 
external loads perpendicular to the xy plane and applied 
moments. [Ref. 5} 

The stresses shown in Figure 6 result in the following 
equations {Ref. 5] for bending moments M and transverse shears 
Q. 


t/2 c/2 e/2 
s = = 10a 
M, f none? My, f ad My, f pp tyeaz (10a) 


0, = f mca Q, = f yada (10b) 

Stresses 6, and oy, are greatest at the surface 2=it/2, 
while +,, is maximum at the midsurface. Transverse shear 
stresses t,,, t,, are small compared too,, 6, andt,,, and are 
not considered in the classical Kirchhoff plate theory. 

In what follows, Kirchhoff’s plate theory [Ref. 5] is 
briefly reviewed, which forms the basis for the GIFTS 
formulation. 

As transverse loads are applied to the plate, the points 
on the midsurface move only in the z direction. Under loads, 
normals to the midsurface are assumed to remain normal before 
and after deformation. Figure 7 shows a differential element 


of a thin plate before and after loading. As shown, the line 


OP is perpendicular to the midsurface before and after loading. 








fe) 


Figure 7. Differential Element of a Thin Plate a) Before 
Loading b) After Loading 


P F Ow ow : 
Bcd = a 
From Figure 7, with ax and ay being assumed as 
small angles, the following relations apply: 
eg lh du &w 
u -2-—= @. = -———— 
ox hence * Ox ax? 
ow av ew 
vz-g = OY =z ez 11 
ay Sy ay * 2 Sa ou) 


- Ou, OY . Lop ew 


Ya = Oy * Ox ~~? Bxay 


These are the strain displacement relations in 





accordance with classical Kirchhoff's plate theory applied to 
a@ thin plate. 

By using stress-strain relations, the moment-curvature 
relations may be derived. Neglecting thermal expansion and 


assuming isotropic material, stress-strain relations are given 


by: 





tv 0 lye, 


iy) = ot vi 0 & (12) 
1-v i-v 
00 =| Wx, 


By substituting equation (11) into equation (12) and 
substituting the result into equation (10a), the following is 


obtained: 
{Mm} = -[D,] {x} (13) 


where the moments and curvatures are 





Sw ew, ew 
= PS = |Z ow 14 
UM aL MM FG ted = LS St toy! (14) 
and 
lv 0 
ina > «2° (15) 
09 -h:v) 
2 
3 
where D= —2& yp is called the flexural rigidity of the 
12 (1-v?) 


plate and is analogous to bending stiffness EI of a beam. 
The next step is to formulate the element stiffness 


matrix using the virtual work equation: 


SWine = fs {$e}? [B) {Se} dv (16) 


- mee cn ne tne omemnnenmaattnn, 


crue napene woniethigtndae ~ tee 








fe}? =1-2S% ~ 2 Fw 57 Sw 
Gx? Gy? oxy 





and where [£] is given by equation (12). é6w,,, and de are the 
internal virtual work and virtual st-ains respectively. 


For plates, using the strain-displacement relations, 





bityne = f 5) 71D, teh da, where {x}? = Ze ne a £23 (17) 


For an element having N nodes, displacements W is interpolated 


as 


wy. Lard 
# = aya O (38) 


where the nodal dearees of freedom are given by 


{a} “lw gee... a = ou 


Equation (18) is differentiated to yield curvatures 
fg} = f(a (19) 
On substituting equation (19) into equation (17) yields 


BWine = (8d}7 LK) {cd} (20) 


where the element stiffness matrix (K] is identified as 


f (2172, Waa (21) 


shen . 





B. FREE VIBRATION OF MULTI-~DEGREE-OF-FREEDOMN SYSTEMS 

Craig (Ref. 6] provides a good overview of muiti-degree-of- 
freedom (EDOF) systexs. 

The equation of motioa for a free undamped MDOF system can 
be written as 

{mj {a} + [K] {u} = {0} (22) 

where [m] and [K] are (NxN) matrices and {u(t)} is a Nxl 
vector o£ generalized displacement cooreinates. The soivtion 


of the differential equation gives harmonic motion Given by 


u =u, cos (we-a) (23) 


Substituting equation (23) into (22) yields the eigenvalue 


problem 


(LX) ~w? [a]) fu} = 0 (24) 


For non-trivial solution, 


Ix - wml = 0 (25) 


Equation 25 is recognized as the characteristic equation for 


the free vibration response. The resulting polynomial in w? 


yields the roots cr the eigenvalues which correspond to the 


natural frequencies of the system. 
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Corresponding to each eigenvalue, we, there will be an 


eigenvector or natural mode u,, where 


u,=4°} r=1,2,...,8 (26) i 


Typically, these modes may be scaled by a process called 
normalization and results in modal vectors called normal 
modes. A mode that has been scaled to have a unique amplitude 
will be designated as $6, and will be dimensionless. The modal 


vector corresponding to w, can be written 
By = co, (27) 


where c, is a scaling constant whose units are such that $,1m $, 
has the dimension of mass. 
There are three procedures for normalizing modes for MDOF i 


systems. [Ref. 6] 


1. Scale the r™ mode so that (6,),=1 at a specified \ 
coordinate i. 

2. Scale the rt mode so that (,),=1, where 
1 @,),1=max,} ,),1, 

3. Scale the r“ mode so that the generalized mass is 
defined by 





cree = ates be neh ntl, 


: n, = 6%mb, (28) 


The generalized stiffr +3 for the x™ mode is 


K, = $2kb, (29) 
By expressing equation (24) for the r' mode and premultipling 
by $ Fae the generalized stiffiness-nass relations are obtained 


as 


(30) 


C. NUMERICAL EVALUATION OF MODES AND FREQUENCIES OF MDOF 
SYSTEMS 


This section discusses the procedure used by GIFTS to 
obtain numerical solution to large eigenproblems. The LACE 
spacecraft in its fully deployed state was modeled by over | 
16,000 degrees of freedom. The dynamic analysis of this 
structure involves determining the natural frequencies and the 


corresponding natural modes by solving the equation 4 


(K-w?) }; =0 (34) 
Vector iteration methods are simple and elegant for obtaining 
eigenpairs. The specific method used by GIFTS is the subspace 
iteration method. The subspace iteration solution is very 
effective in the calculation of the lowest eigenvalue and 
corresponding eigenvectors of systems with large bandwidth and 


which are too large for the high-speed storage of the 


computer. (Ref. 7] 
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Bathe [Ref. 8] establishes three steps for the subspace 
iteration method. 


1. Establish q starting iteration vectors, qg>p, where p 


is the number of eigenvalues and vectors to be 
calculated. 


2. Use simultaneous inverse iteration on the q vectors 
and Ritz analysis to extract the best eigenvalue and 
eigenvector approximations from the q iteration 
vectors. 

3. After convergence, the Sturm sequence check is used 
to verify that the required eigenvalues and 
corresponding eigenvectors have been calculated. 

The objective of the subspace iteration method is to solve 


for the lowest p eigenvalues and eigenvectors satisfying 
Ko = MOA (32) 
where A =diagonal (A;) and $={$),-.- -6,]. The eigenvectors must 
also satisfy the orthogonality conditions 
OTK = A; OTM =F (33) 
Detailed derivationg of the subspace iteration method are 
shown in Bathe (Ref. 8]. The subspace iteration algorithm 
shown below finds an orthogonal basis of vectors in E;,, 
subspace. 
Xin = MX, (34) 
for L=1,2,..., and with iterations from E, to E,,,. Next, 


projections of the operators K and M onto K,,, are computed: 
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Ky. = Kea KX9 (35) 


nate 
. 
SAF SRN Sa enema, 


My = XoMXy 135) 


By solving for the eigensystem of -projécted. operators 


SA ASee 8 
we 


Kis Quer 2H, DriaAa (37) ; 


where Q is an orthogonal matrix, improved approximation for . 
the eigenvecters is found -by 

Sees = sae Or01 : (38) 
It may be aoted that A,,,--> A ang X2, ~~? $ aS L-~> &, 

The first step of the suhspace iteration is to generate the 
starting iteration vectors in %,. The following algorithw is 
used to select the starting iteration vector. The first column 
in MX, is the diagonal of M. The other columns are unit 
vectors with entries +i at coordinates with the smallest k,/m, 
ratio. 

The subspace iteration method requires a measure to compute 


convergence. Assuming that in (L-1) and L iterations, 


eigenvalue approximations 4,4 anda "+), is1,...p, have been 


calculated. The measure for convergence, then, is 


pases? - a fay 


- f= 3 
ry stol ; i=1,...p. (39) 


where tol may be 10°5, when eigenvalues shall be accurate to 


2S digits. 





Since equations (32) and (33) can be satisfied by any 
eigenpairs, there must be a way to verify the calculations. 
Once the convergence is satisfied in equation (39), with s 
being at least equal to 3, a check may be performed to make 
sure that the smailest eigenvalues and corresponding 
eigenvectors have been calculated. The Sturm sequence property 
is used to provide this check. This property is derived from 
the following analysis. [Ref. 68} By using the Gauss 
elimination solution, the stiffness matrix can be factorized 
as 

K=LDL? (40) 
where L is a lower unit triangular matrix and D is the 
diagonal matrix. 

Let K-+1M be factorized into LDL’. In the decomposition of 
K-1M, the number of negative elements in D is equal to the 
number of eigenvalues smaller than pp. Because of this 
property, by assuming a shift p and checking whether p is 
smaller or larger than the required eigenvalue, successive 
iterations can reduce the interval in which the eigenvalue 


must be. A summary of subspace iteration solution is shown in 


Table II. 
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Table II. SUMMARY OF SUBSPACE ITERATION 


III. LACE FINITE ELEMENT MODEL 


A. GIFTS CAPABILITIES 

GIFTS has several capabilities which facilitate the 
modeling of large complex structures. Some of these 
capabilities include: 

1. automatic model generation 

2. model editing, display and information 

3. automatic load and boundary conditions 

4. vibrational mode extraction 

5. substructuring 

6. thermal stress analysis 

Element types that can be used include: 

1. rods and beams 

2. plates/shells 

3. solid elements and axisymmetric elements [Ref. 9] 

The elements can be selected from a library of options. The 
materials can be created by the user or from a library of 
defined materials. GIFTS also allows users to define 
anisotropic materials. 

Static and dynamic analysis can be performed on the model. 
The dynamic analysis provides free vibrations and mode shapes. 
These vibrational mode shapes can be displayed on the screen. 


For structures that undergo thermal loading, deflections and 
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stresses can be calculated by computing thermal (pseudo- 
forces) forces to simulate the thermal effects. 

GIFTS also has the capability for substructuring and multi- 
level substructuring. Large models may be divided into 
substructures. Certain areas of a substructure may be modeled 
as a second level substructure to allow economic modeling and 


reduce computational costs. 


B. SIMPLE FINITE ELEMENT MODEL OF LACE 

Initially, the LACE spacecraft is modeled as a point mass 
with three attached beam members. The three dimensional 
triangular trusses are modeled as solid circular beams. To 
ensure that this beam had the same bending characteristics, 
the following bending stiffness relations was used from A=C- 


Able Engineering (Ref. 10). 


flexural rigidity (EI) =1.5nER‘e (41) 
where 
e = maximum bending strain of longerons when completely 
coiled @ = d/2R = F/E) 
F = coiling stress of longerons 
ad = longeron diameter 
E = Young’s modulus of longeron material 
R = boom radius 
with 
e = .015 
R = 5.0 inches 
E = 8.0 x 10° psi 


This results in EI = 5.3 x 10° lbs-in’. 
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From this relation, I is seen to be .663 which results in 
a solid circular beam with a radius of .9583 inches. The beam 
is with 10 elements for a 150 foot truss and 5 elements for a 
75 £t. truss. This model is used in GIFTS with the appropriate 
tip masses to determine the natural frequencies and mode 
shapes. The resuits will be shown and discussed in the next 
chapter. A listing of the file to generate the model is given 


in Appendix A. 


C. COMPLEX FINITE ELEMENT MODEL OF LACE 
1. Main Spacecraft Body 

The on~orbit configuration of the LACE spacecraft is 
shown in Figure 8. From the figure, it can be seen that the 
LACE spacecraft consists of structural panels, solar panels, 
sensor panels, truss elements and various beam types. 

The primary structure of LACE consists of a frame type 
structure consisting of channel section stringers, tee-section 
stringers, z-section doublers, and angle-section members. The 
primary structure also consists of honeycomb panels. Figures 
9a and 9b show the basic configuration of the primary 
structure. 

The secondary structure of the LACE spacecraft consists 
of the fixed and deployable solar array substrate, deployable 


sensor panels and deployable sensor arms. Figure 10 shows the 


LACE spacecraft secondary structure. 
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Figure 8. LACE Spacecraft On-Orbit Configuration 


The primary frame members are made of aluminum, AL6061- 
76. Table IIZ shows the material characteristics of AL6061-T6 
used in GIFTS. ' 
Table III. AL6061-T6 CHARACTERISTICS 


2.5382 E-4 lbs-sec*/in’ 


The honeycomb panels consist of .05 thick face sheets 










consisting of aluminum, AL6061-T6. The core of the panels are j 


1/4 -5052 -.003, with a density of .01552 lbs~sec*/in’, 
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Figure 9a. LACE Spacecraft Primary Structure 


Figures ila and 1lb show the primary and secondary 
honeycomb panel structure and dimensions. 

The primary panels are one inch thick with a .9 inch 
core. The solar array panels are .5 inches thick with a .4 
inch core and the deployable sensor panels are .75 inches 
thick with a core of .65. The deployable sensor arms are 
channel beams with .125 inch thick aluminum, AL6061-T6. 

The honeycomb panels are modeled as aluminum panels, 


AL6061-T6, with appropriate conversions to account for 
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Figure 9b. LACE Spacecraft Available Panel Area 


equivalent stiffness and weight. The panel stiffness for a 


honeycomb panel is given by {Ref. 11] 


= —Fth? (42) 
2 (1 = v?) 
where 
E = Young’s modulus 
t = face thickness of the panel 
h = core thickness 
v = Poisson's ratio 


The stiffness for each honeycomb panel is calculated 


using the above formula. 
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Figure 10. LACE Spacecraft Secondary Structure 


t stiffness, an equivalent thickness is 


Using the 
calculated, using the flexural rigidity formula for isotropic 


materials. 
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Figure lla. Primary Panel 
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Figure lib. Secondary Panels 


2 —2t?__ (43) 
12 (1-v?) 


The mass of the honeycomb panels is computed and used 
with the new thickness to determine volume and density. These 


calculations are presented in Appendix B. Table Iv summarizes 
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the new parameters used to modéi thé honeycomb panels as 


aluminum panels with appropriate stiffness and mass. 


Febie IV. THICKNESS AND DENSITY PARAMETERS 


Solar Panels 


Sensor Panels 6.22 x 10° 





The £inite elements model of the main body of the LACE 
spacecraft is shown in Figure 12. The associated file to 
generate the model is given in Appendix C. 

The LACE spacecraft has several sensors and components. 
Even though the spacecraft is essentially a rigid body, the 
mass of the components was modeled as accurately as possible 
according to the mesh size of the grids. Appendix D contains 
the component placements and mass distributions. 

2. Spacecraft Trusses 

The automatic deployable lattice booms are manufactured 
by AEC-Able Engineering Company, Inc. They are designed for 
applications that require high dimensional stability and high 
ratio of bending stiffness to weight. Figure 13 shows the 
principal parts of the continuous longeron boom and the 
retraction geometry. The longerons are continuous along the 
length and are connected to the batten frames with pivot 
fittings. Six diagonals provide shearing strength and 
stiffnesses. When the boom is twisted, tension increases on 


three of the diagonals, causing batten members to buckle and 
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Figure 12. Finite Element Model of LACE Spacecraft 


shorten. As twisting increases, the longerons rotate about the 
pivots and assume a helical configuration. In the retracted 
position the longerons are coiled in flat helices while 
battens lie on top of each other. 

The longerons and battens are modeled as circular beams 


while the diagonals are modeled as rod elements. The battens 








TRANSINONAL 
CONSIGUZ ATION 


Figure 13. Continuous Longeron Boom 


are oval-shaped and attached to the longeron by pivot joints, 
but for simplification the battens are modeled as circular and 
the joints are not modeled. Listed in Table V are the 
dimensions and properties of the truss. When fully retracted 
the triangular batten frames lie within a 10 inch diameter 
circle. 

Figure 14 shows a cross-sectional view of the truss and 


its dimensions. 
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Table ¥V. TRUSS PROPERTIES AND DIMENSIONS 


Longeron and Batten Diameter -150 inches 


Diagonal Diameter -050 inches 
Material Type S-glass epoxy 
Density -075 lbs/in? 


Young’s Modulus E=8.0 x 10° psi 





Figure 14. Cross-Sectional View and Dimensions of Longeron 
Structure 


The booms of the spacecraft are ideal for 
substructuring. The 150-foot boom results in over 5,000 
unknowns, when modeled as three-dimensional beam-truss 
elements. The suhstructuring technique is used in modeling 
when the number of unknowns may be reduced substantially. 


Initially, substructuring and multi-level substructuring 
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techniques were used to model the trusses, reducing the 
computational time considerably. 

The first substructure was a 15-foot section of the 
truss. Ten of these substructures were joined together to 
construct the 150-foot boom. This multi-level substructure was 
attached to the main body of the spacecraft, thus forming a 
super element to the main body. This reduced the complete 
model from over 16,000 unknowns to about 2,500 unknowns. This 
was an attempt to capture the dynamics of the LACE spacecraft 
more accurately. 

However, the use of substructuring generated large 
negative masses in the mass matrix. Kamel et al. (Ref. 12] 
provide a detailed formulation of the constrained 
substructurang techniques. As only the executable version of 
GIFTS program was available, and the objective was to obtain 
dynamic characteristics, it was decided to pursue modeling the 
whole structure. 

The limitations of the program or the methodology 
adopted there is being pursued as well. The supplier of 
software is looking into the problem. The input files used to 
generate the multi-level substructuring model of LACE and 
typical negative mass elements are attached in Appendix E. 


The listing of the file to create the trusses and 


attachment to the spacecraft is in Appendix F. 





Iv. RESULTS OF DYNAMIC ANALYSIS 


For linear behavior, resonance occurs when the frequency i : 
of excitation equals the natural frequency. In order to avoid 
the ill-effects of large amplitude vibration at resonance, the 
natural frequency must be known and compared with potential 
excitation frequencies. The gravity gradient pitch vibration 
frequency is 2.3 x 10% Hz and is widely separated from the 
lowest modes of the finite element models. ; 

Table V presents computed values for the first four modes 
of three different models of the LACE spacecraft. 


Table V NATURAL FREQUENCIES OF THREE LACE MODELS 


Mode GIFTS NASTRAN [Ref 1] GIFTS | 
Beam Model ae Beam Model ae ae Mode} (Hz) 


1.01930] | 01938 | faa | 


‘ ' 
= 04825 £04729 .0516 


-05454 -0536 -0588 


ma 
es re 


Initially, the NASTRAN beam model was developed by Naval 





Research Laboratories (NRL). The present beam model was 


developed as described in the previous chapter. The first 


' 
| 
three modes agree within 2%, however the present fourth mode 
appears to be an anomaly. The GIFTS complex model is 


consistently 10% higher than the NASTRAN model. The NASTRAN 


model appears to yield fairly good data in the lower modes. 





This discrepancy may be attributed to the modeling 
uncertainties and approximation of the geometric and stiffness 
distributions. The detailed modeling of NASTRAN was not 
available. It may be noted that there is no excitation 
frequencies at those computed frequencies. Further, complex 
modeling is recommended when higher modes and frequencies are 
required. The simple model may not capture these higher modes 
and even skip some modes. The higher modes assume importance, 
especially, for very flexible structures and in the design of 
control systems for vibration control and suppression. Figures 
15 to 20 show the mode shapes and frequencies for the NASTRAN 
model while Figures 21 to 26 show the data for the GIFTS beam 
model. Figures 27 to 30 show the GIFTS complex model of the 
spacecraft. Appendix F contains the frequency and mode shapes 


of the LACE spacecraft in three configurations as it deploys 


to its final state. 
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-04729 Hz, NASTRAN BEAM MODEL 


Figure 16. Mode 2, u 








-05362 Hz, NASTRAN BEAM MODEL 
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Figure 17. Mode 3, v 


-1104 Hz, NASTRAN BEAM MODEL 


Figure 18. Mode 4, vu 
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-1808 Hz, NASTRAN BEAM MODEL 








Figure 19. Mode 5, u 


20s 00 on ee em on oot on em om om ee we 
. 


optic a a es ae tes 





NASTRAN BEAM MODEL 
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Mode 6, v = 


Figure 20, 
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Figure 21. Mode 1, v = .01930 Hz, GIFTS BEAM MODEL 








Figure 23. Mode 3, v = .05459E Hz, GiFTS BEAM MODEL 





i 
JA 
4 ¥ 
t 
. 5 
' ' 





Figure 24. Mode 4, vu = .1738 Hz, GIFTS BEAM MODEL 
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Figure 25. Mode 5, v = .2142E Hz, GIFTS BEAM MODEL 
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Figure 26. Mode 6, v = .2349 Hz, GIFTS BEAM MODEL , 








Figure 27. Mode 1, u = .02163 Hz, GIFTS COMPLEX MODEL 








Figure 28. Mode 2, v = .05165 Hz, GIFTS COMPLEX MODEL 








Figure 29. Mode 3, v = .05879 Hz, GIFTS COMPLEX MODEL 
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V. THERMOELASTIC EFFECTS 


A change in temperature along a bar will change its 
dimensions. When an isotropic bar is heated uniformly and is 
free to expand, the sides will increase in length. The 


material undergoes a uniform thermal strain e, given by: 
e, = a(AT) (44) 


where @ is the coefficient of thermal expansion and AT is an 
increase in temperature. The length of the bar will increase 


by an amount 
8. = a(AT)L 


where L is the length of the bar. 

The ABLE booms (Ref. 10} used on LACE are designed so that 
they undergo minimum thermal bending or twisting in the solar 
radiation environment. Pretwist is used to prevent thermal 
twisting or thermal Lending that would occur if one longeron 
18 shadowed by another. 

This chapter presents analysis for deformations that could 
result in a worst case scenario. Two possible scenarios 
considered are when the boom may bend due to unequal heating 
of the diagonals and unequal heating of the longerons. Unequal 


heating of the diagonals is more likely to occur than unequal 


heating of longerons. [Ref. 13}. The following analysis 
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cavetiene only the unequal heating of the diagonals and 
assumes no shadowing effects. 

Sources of heat include spacecraft components, solar flux, 
albedo flux and thermal radiation of the earth. The solar flux 
is defined as the flux existing at a distance of one 
astronomical unit (AU) from the sun. Albedo flux is the 
fraction of total incident solar radiation on the earth which 
is reflected into space as a result of scattering in the 
atmosphere and reflection from the clouds and earth surfaces. 
Thermal radiation from the earth is the portion of incident 
solar radiation absorbed by earth and its atmosphere and re- 
emitted as thermal radiation according to Stefan-Boltzman law 
(Ref. 11]. For the computation of LACE thermal deformations, 
the following data is used {Ref. 13}: 

Solar flux: 442 Btu/sq ft hr 

Thermal Garth Radiation: 70 Btu/sq ft hr 

Albedo flux: 160 Btu/sq ft hr 

e/a = .8 
where e is the emissivity anda is the absorptivity. 

The worst case of unequal diagonal heating occurs when the 
sun rays are parallel with one set of diagonals and almost 
perpendicular to the other set. The set of diagonals 
perpendicular to the sun will receive maximum solar flux. The 
parallel set will receive no solar flux, but will receive 
earth’s albedo and thermal radiation flux. A simple approach 


is taken where the hottest and coldest temperatures cf the 
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diagonals are calculated. This results in a maximum 
temperature of 192°F and a minimum of -5°F. Following the 
discussion in Ref. 12, calculations are carried out and 
presented in Appendix G. 

It is assumed that the battens and longerons receive equal 
heating of 80°F. Using a 3-D model of the truss, the 
temperatures above, and a coefficient of thermal expansion of 
1.75 x 10°°/°R, the maximum deflection is calculated. Figure 
3la shows that the deformed shape is similar to the first 
bending mode. Figure 31b shows a closeup of the deformation. 
A maximum deflection of 1.88 inches is calculated using GIFTS. 
This is a 1.2% deformation for the 150 foot truss and should 
have negligible effect on the system dynamics. 

However, thermal flutter could occur if the period of the 
thermal loading coincides with the natural frequencies of the 
spacecraft. 

The first bending mode has a natural frequency of .02163 
Hz. The thermal loading from the sun occurs over a long period 
and is on the order of 10%. The thermal loading is two orders 


of magnitude lower than the fundamental frequency and as a 


result should not interact with any natural frequencies. 
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Figure 3la. poeta elastic effects resemble 1st bending mode 
shape. 
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Figure 31b. Close-up view of thermoelastic effects 





VI. MULTI-BODY DYNAMICS 


A. COMPCNENT MODE SINTHESIS 

In the previous chapter, finite element techniagues were 
used to formulate a model of the LACE spacecraft for use in 
structural dynamics analysis. This section discusses a class 
of reduction methods known as component mode synthesis, or 
substructure Coupling for dynamic analysis. These methods are 
usezul for analysis of large structural dynamics problems. 

The basic idea of component mode synthesis is to treat the 
complex structure as an assemblage of substructures. Each 
substructure is analyzed independently and then their dynamic 
characteristics (mode shapes and natural frequencies) are 
synthesized to analyze the complete structure. There are many 
variations of the method of component mode synthesis and 
extensive iitera-ure is available (Ref. 14, 15, 16, 17). 

Hurty (Ref. 14] developed a procedure for analysis of 
structural systems using a displacement method which used 
three types of generalized coordinates, namely: 1) rigid body 
coordinates, 2) constraint coordinates, and 3) normal mode 
coordinates. Hucty used Rayleigh-Ritz approach in his 
formulation. The Craig-Bampton method [Ref. 15] is similar to 


the treatment due to Hurty, except that it simplifies the 


treatment of rigid-body modes of substructures by eliminating 





the separation of boundary forces into statically determinate 
and statically indeterminate reactions [Ref. 15]. MacNeal 
(Ref. 16] introduced residual flexibility modes to retain the 
static contribution of a higher frequency truncated modes. 
Rubin [Ref. 17} developed a new method which adds residual 
inertial and dissipative results to the method introduced by 
MacNeal [Ref. 17}. 

This section will present the basics of component mode 
synthesis. 

Craig (Ref. 18] provides a good overview of component mode 
synthesis methods. His notation and examples will be used 
extensively in this discussion. 

A substructure is generally connected to one or more 
adjacent components and is composed of interior degrees of 
freedom and boundary degrees of freedom. Figure 32 illustrates 
a substructure connected to other components and shows 
boundary and interior coordinates. 


The equation of motion for a component is given by 


MK + CX + KX =f (45) 
where M is the mass matrix, C is the damping matrix and K is 
the stiffness matrix. 

The total set of physical coordinates of the component is 
defined as P, while interior and boundary points will be 
defined as I and B respectively. Boundary coordinates are 


further subdivided into statically determinate (rigid body) 


62 





te 


e Boundary o Interlor 
coordinate coordinate 





Figure 32. Unrestrained truss component 


support of the component and its complement, redundant 
(excess) boundary coordinates designated as R and E 
respectively. These designations are used throughout this 
analysis. 

Figure 33 shows a beam divided into several components and 
will be used to to illustrate mode sets in the following 
discussion. 

The two fundamental steps of component mode synthesis are 
to: 1) define component modes, and 2) to define coupling of 
components to form a system. 

The physical coordiaates, x, can be represented by 


component generalized coordinates, P by the transformation 


X= yP (46) 
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Figure 33. Beam model divided into components 


where ¥ consists of component modes of specified type [Ref. 
14]. These include: 

normal modes of free vibration 

attachment modes 

constraint modes 

rigid body modes 

inertia relief modes 
These modes are defined as follows. [Ref. 17) 

1. Normal Modes 

Normal modes are classified as fixed interface normal 

modes, free-interface normal modes or hybrid-interface normal 
modes depending on how the interface coordinates are 


restrained when the component normal modes are obtained using 





Cee any 


seem me acinomae 


= team 


(K - Wgm)op = 0 (47) 


The modes are normalized with respect to the mass matrix m. 


$209, = Ing $2K, = Aj, = diag(Wp)? 
where 6, is component normal modes. The notation 6 will be 
used for normal modes, while y will be used for assumed modes. 
2. Constraint Modes 
A constraint mode is defined by statically imposing a 
unit displacement on one coordinate of a C set of physical 
coordinates. Let C=E and define a constraint mode by placing 
a unit displacement on one coordinate of the C set and zero 
displacement on the remaining C set. The matrix of constraint 


nodes, ¥, is defined by the equation: 


Kez Kge Kiel [Wye ° 
Ke Kee Ky Log = Ree (48) 


Key Kr K,, ° Re! 


This equation may be simplified to yield 


-1 
Pie 
“Kj, Ky 
Ye, = Zee = I 7 (49) 
fe) ce 
e) 


3. Attachment Modes 
Attachment mode is defined by applying a unit force of 
the coordinates of an A set (Ref. 15}. In this case, 
attachment modes will be defined for A = E. The matrix of 


attachment modes », is shown by the equation 
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Kiy Kia Kiz| [¥,.] [0 
Kaz Kaa Kaz Weal = Taz (99) 
Ki; Kya Kr ° Rey 


Using the first two rows of equation (50), %, may be 


represented by 


Wie Gia 
Wa = [Vaal = |9aa (51) 
oO fo) 


where g = K! is the flexibility matrix. 
4. Rigid Body Modes 
The boundary conditions are depicted in Figure 45, where 
the R set will restrain the component from rigid body motion 
and the E set contains redundant boundary condations. By 
defining the’ ragid body modes relative to the R set, the 


equation is given by 


Kis Kie Kye] [War [o) 
Kir Kee Ker Wor = 10 (52) 
Key Kre Key Ty, © 

which simplifies to 
Ky; a Mg < fe (53) 
Ky: Kee. er, kK, 





The y, matrix is given by 
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(54) 
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‘i 5. Inertia-Relief Modes 
This method defines attachment modes for a component 
with rigid body freedoms. This method was presented by Rubin 
{Ref. 16] and MacNeal (Ref. 17}. By letting D’Alembert force 
vectors associated with rigid-body modes be applied statically 
to a component which is fully constrained on the boundary, $, 


can be defined as follows. 


Kis Xto Kizl Wy) [Mat Mee Miz} War} | 0 
Kos Kee Ker} | O | = |Moxr Mos Merl [Wor] + rer x25) 
Ky K,, Kr; ie) My Mee Mr, Ty Re 


Var] [KO (My War = MieWor = Mir) ; 
%, =] 0] = fe) (56) 
fe] (e] 


6. Coupling of Components 
This section describes generalized substructure coupling 
as applied to free vibration analysis. [Ref. 18] 
Assuming two componentsa andf having a common boundary 


interface, compatibility of interface displacement requires 


xf = xt (57) 


The interface forces are related by 
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ff = fh =0 (58) 
By representing the physical coordinates x, by 


generalized coordinates p, the following equations are derived . ‘ 
xi = pepe, x8 = whee (59) i 


where y* and y® contain assumed static and dynamic modes. 
The constraint equations can be written in generalized 


form to form a single constraint equation 


CP=0 (60) 


: 61 
oP (61) 


Let P be rearranged and partitioned into dependent P, and 


where 


independent, P,, coordinates. Then, 


P, 
{Cop Coz) fe 20 (62) 
Py 


where C,, is nonsingular square matrix and equation 


P, = Sg (63) 





defines S and g, where 


wor 
2 -| Cop =| (64) 
Try 


The p andx corresponding to P are given by 





a x 

pet ol. x =ft ‘4 (65) 
oO pel o xb 

fhe coupled system of equations for an undamped system 


is given by 
Mg + Kg =0 (66) 


where 


M=S™S, K= SS (67) 


B. DYNAMICS OF FLEXIBLE BODIES IN TREE TOPOLOGY 

1. Overview of Multibody Systems 

Spacecraft and large spacecraft structures are typical 

multibody systems. Large strides have been made in the last 20 
years in the efficient formulation and solution of multibody 
systems. Particular interest in multibody dynamics has risen 
in spacecraft dynamics. Initially, space vehicles were 
idealized as rigid bodies or elastic beams. In the mid-1980s 
equations of motion were published for a point-connected set 
of interconnected rigid bodies in a topological tree. A model 
containing rigid bodies and elastic appendages was eveloped in 
the 1970s. The next major step was the incorporation of body 
flexibility in the topological tree model. (Ref. 19] 

2. Multibody Computer Program - TREETOPS 


TREETOPS is a computer program developed to deal with 


multibody structures in an open-tree topology. It is a time 





history simulation of a complex multibody flexible structure 
with active control elements. Some of the features include: 1) 
any or all bodies can be rigid or deformable, 2) hinges can 
have zero to six degrees of freedom, 3) the dimension of the ; i 


problem equals the number of degrees of freedom, 4) individual 


ne 
1 
eb ea reicana e 
w= tenet ememaennenesaneanneanmmtans, 
+ 
Re ee a arr etnnemernnn cena 


body deformation can be described by any set of modal vectors, 
5) an interactive program, and 6) extensive control simulation 
capability. (Ref. 18) 

The computer simulation consists of three parts: 
defining a tree topology of flexible structures, define a 
controller and a set of sensors and actuators. 

The structure of the TREETOPS model consists of bodies 
and hinges with sensors and actuators included for interfacing 
with the control system. Figure 34 shows a typical structure. 

a. Body Types 
The program simulates a set of bodies in a tree 
topology. Each body is defined independently. Sensors, 


actuators and hinges are connetted to specific points called 





node points. Each body may be defined as rigid or flexible. 
Each body is to be defined with an ID number, mass 
properties, center of mass and all its pertinent node points. 
(Ref. 20) 
b. Hinges 
Hinges interconnect two adjacent bodies. One body 
is called the inboard body and the other is called the 


outboard body. The inboard body is the one on the side closest 
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Figure 34. Structure composed of bodies and hinges 


to the inertial reference frame. Figure 35 shows a hinge 


representation. 


The functions of the TREETOP hinges are to 
(Ref.20): 


1. define topology of the structure 
2. define kinematic variables of the multibody system 


3. define relative orientation between adjacent bodies. 


c. Sensors and Actuators 
A set of 16 sensors have been buil< into the 


simulation, They include rate gyros, resolvers, 
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Figure 35. General representation of jth hinge 


accelerometers, position and velocity sensors, tachometer, sun 
and star sensors, etc. The actuators serve as a way to apply 
force and torque inputs. Inputs may be control or disturbance 
inputs. Disturbance inputs can be applied with function 
generators to the actuators. Actuator types include reaction 
jets, hydraulic cylinder, moment actuator and torque motor. 
d. Orbit Environment 

TREETOPS has the capability to model the orbit 
environment of a spacecraft to include gravity gradient and 
aerodynamic drag. A magnetic field model is included which 
produces a force through interaction with magnetic actuators. 


In computing the atmospheric drag, TREETOPS uses atmospheric 








density comouted from the Gacc 


spacecrs it orbit is Get 


parameters. 


HIB 3970 atrospheric model. 


. 
=e 


cmined by entering the six orbital 
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VII. CHAOTIC VIBRATIONS 


A. HOW TO IDENTIFY CHAOTIC VIBRATIONS 

Chaos is defined as a motion that is sensitive to initial 
conGitions [Ret. 21]. Chaos can occur oniy in nonlinear 
systens, but all nonlinear systems do not exhibit chaos. Chaos 
can be observed in many physical systems. A partial list is 


shown: 


. Vibrations of buckled elastic structures 

2. Mechanical systems with play 

3. Large, three-dimensional vibrations of structures 

4. Kereoelastic problens 

5. Systems with sliding friction 

6. Feecback control devices 

In order to identify chaotic motions, several procedures 
are suggested (Ref. 21], such as 

1. Identity nonlinear elements in the system 

2. Check for sources of random input in the systen 

3. Observe time history of the measured signal 

4. Observe phase plane history 

5. Examine Fourier spectrum of the signal 

6. Determine Poincaré map of the signal 


7. Vary system parameters 
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Quantitative characteristics of chaotic vibrations and the 
diagnostic tools used are summarized as follows: 


1. Sensitivity to changes in initial conditions 
(Lyapuncv exponent and fractal basin boundaries) 


2. Broad band spectrum of Fourier transform when rxotion 
is generated by a single frequency 


3. Fractal properties of the motion in phase space which 
indicate a strange attractor (Poincaré maps, fractal 
dimension) 


4. Increasing complexity of regular motions as 
experimental parameters are changed 


5. Transient or intermittent chaotic motions; 
nonperiodic bursts of irregular motion or initialiy 
yandom-like motion that settles into requiar motion. 

1. Nonlinear System Elements 
A linear system does not exhibit chaotic vibrations. 

fypical nonlinear effects from mechanical systems include 
nonlinear stiffness, material nonlinearity, nonlinear damping, 
free-play, and nonlinear boundary conditions. Nonlinear 
elastic effects can be due to large deformation. A good 
example of material nonlinearity is the stress~strain 
relations of materials modeling rubber or elastomers. 

2. Random Inputs 
There are no assumed random inputs in chaotic 

vibrations. Applied forces and excitations are assumed to be 
deterministic. By definition, chaotic vibrations arise from 
deterministic physical systems. A large output signal to input 
noise ratio is required if nonperiodic response is to be 


attributed to a deterministic system behavior. [{Ref. 21) 
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3. Observation of Time History 

The first indication of chaos may be indicated in the 
time history. The motion observed shows no visible pattern or 
periodicity and may be chaotic or random. 

This method is not conclusive, since motion could have 
a long-period behavior that is not detected or quasiperiodic 
motion where two or more periodic signals are present. 

4. Fourier Spectrum 

The presence of a broad band Fourier spectrum in the 
output is another clue that may be used to suspect chaotic 
vibrations. A precursor to chaos is the presence of w,/n 
subharmonics. However, multiharmonic outputs do not always 
imply chaotic vibrations as hidden degrees of freedom may be 
present. 

5. Phase Plane History 

in the phase plane, complete information about a 
Gynamical system is represented by a point. At the next point 
in time when the system dynamics change, the point is 
displaced. 

This moving point gives the history of the dynamical 
system. The coordinates chosen for the study of dynamics, 
typically, are the amplitude and velocity of motion. Figure 36 
shows the phase plane of a simple pendulum. The circle on the 
phase plane represents the motion over one cycle and is called 
the trajectory. It should be noted that a periodic motion is 


a closed orbit in the phase plane and is called a limit cycle. 
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Figure 36. Phase plane of a pendulum [Ref. 21] 


Chaotic motions have orbits that never close or repeat. As a 
result, the trajectory of the orbits will tend to fill up the 


phase plane. 
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6. Pseudo-Phase-Space Method 

This method is used when only one variable is available, 
as in typical flicht tests where measurements are from a 
strain-gage or accelerometer. In order to do the 2-D phase 
plot from strain-gage measurements, the signal must be 
differentiated. In the case of accelerometer data, the signal 
has to be integrated twice. 

However, by integrating or differentiating, the signal 
is filtered (Ref. 22]. Differentiating the signal will amplify 
high frequencies and attentuate low frequencies. Integration 
will have an opposite effect. As a result, phase plane plots 
obtained from experimental data will be inaccurate. This 
resulted in the development of the pseudo-phase-space method 
or embedding space method. For a one degree system with a 
measurement x(t), the signal is plotted against itself but 
delayed or advanced by a fixed time constant: (x(t), x(ttT)]}. 
This plot yields properties similar to the classical phase 
plane. The closed trajectory in the classical phase plot wili 
be closed in the pseudo-phase method and chaotic motion appear 
chaotic in both phase planes. 

When the state variables are greater than three 
(position, velocity, time), higher dimension pseudo-phase- 
space may be constructed using multiple time delays, i.e., 
(X(t), x(t+P), x(t+2T)). 

The advantage of the pseudo-phase plane method is that 


a single observable variable can be used to construct the 
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pseudo phase picture and portray the system dynamics without 
distorting the response through integration or 
@ifferentiation. 

7. Poincaré Section 

The Poincaré section can be constructed by placing a 
two-dimensional surface in a three-dimensional phase space and 
noting where the points of the trajectory penetrate this 
surface. This slice will reveal the internal structure of this 
location. 

If the Poincaré section does not consist of a finite set 
o£ points or a closed orbit, the motion may be chaotic. For 
some lightly damped systems, the Poincaré section of chaotic 
motion appears as a set of unorganized points. This motion is 
called stochastic and is shown in Figure 37a. In damped 
systems, the Poincaré section appears more organized with 
parallel lines as shown in Figures 37b and 37c. The Poincaré 
sections can be enlarged to observe further structure (Figure 
38). After several enlargements, if the structured sets 
continue to exist, the motion is defined as a strange 
attractor. This embedding of structure within the structure 
indicates fractal nature of the behavior, which is a strong 


indicator of chaotic motions. 


B. QUANTITATIVE TESTS FOR CHAOS 
The previous section summarized qualitative methods that 


require experience to evaluate chaotic systems. There are some 
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Figure 37. Poincaré maps chaotic motion 


quantitative measures to study chaotic motions. Two well-known 
methods are the Lyapunov exponent and fractal dimension (Refs. 
22, 25]. These methods are described below. The Lyapunov 
exponent will not be used in this analysis and will be 
discussed summarily. 
1. Lyapunov Exponent 

The Lyapunov exponent measures how sensitive the system 

is to changes in initial conditions. It measures the 


exponential attraction or separation, of two adjacent 
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Figure 38. Poincaré map of chaotic vibration 


trajectories in phase space with different initial conditions. 


It is defined as \ 
d(t) = dt (68) 
or 


L= log, (SA) (69) 


where dis the initial distance between two trajectories. 
dct) is the distance at a later time. 


L is the Lyapunov exponent. 
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A positive exponent implies d(t), the later distance 
will be larger than the initial distance and indicates chaotic 
dynamics. 

2. Fractal Dimension 

fhe fractal dimension is another quantitative test for 
chaos and gives a lower bound on the number of essential 
variables needed to model the dynamics of the system. Non- 
integer values for a fractal dimension indicates presence of 
a strange attractor. [{Refs. 22,25] 

There are six ways to classify fractal dimensions. The 
dimension that will be discussed in this analysis is the 
correlation fractal dimension. 


The correlation fractal dimension is defined as 


C(r) =r4 (70) 
where: C(r) is the probability of the attractor within a 
circle, sphere or hypersphere of radius r, and d is the 
fractal dimension. 
By taking the natural logarithm of both sides of 


equation (70) and solving for d, following equation results: 


.,.. {in ¢(r) 11 
ed enc } ee 


The procedure adapted in Sarigul-Klijn [Ref. 22] is 


described below: 


1. Start with a point on the attractor and calculate the 
number of points inside a circle of radius r. 





2. Calculate probability C(r) by dividing this number of 
points by total number of points in the attractor. 


3. Repeat this for several points along the attractor. 
4, Compute C(r) for several values of r. 


5. The slope of log (C(r)) versus log (r) gives d, the 
correlation fractal dimension. 


To obtain the correlation dimension of the attractor of 
a given system, the procedure must be applied in the pseudo 
phase space for several embedded dimensions. The asymptotic 
value of the correlation dimension is the fractal dimension of 


the attractor and is given by 


NON 
1 
Clr) = lim — H(r - |x, - x5|) (72) 
lin op LY bes 351 


where: H(s) = 1 if s>0 and H(s) = 0 if s<0. 

|x, - x,| is the Euclidean distance between the points. 

N is total number of points. 

if the fractal dimension is approximately equal to the 
phase space used for the calculation, the attractor lies ina 
higher dimensional phase space. If the fractal dimension is 
non-integer and is independent of the dimension of phase 


space, the signal is characterized as chaotic. (Ref. 22] 


C. LACE SPACECRAFT BOOM AS A NONLINEAR SYSTEM 
The booms for the LACS spececraft have a constant EI 
distribution [Ref. i0] and the response would appear to be 


linear. However, at each bay, battens are joined to the 
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longerons by joints. Freeplay introduced at the joints tend to 
intreduce nonlinearities. In this section, the spacecraft boom 
élone is modeled and a parametric study to simulate the 
influence of nonlinearities introduced by the freeplay at the 
joints on the system response is presented. The resulting 
behavior of the motion is studied using the methods of chaos 
(Ref. 21}. 

Tae LACE spacecraft boom is modeled as a nonlinear single- 
degree of freedom system. The equation of motion for such a 


system is described by 
mu+f(u, u, t) = A, coswt (73) 
where £(u,4,t) contains nonlinear damping and stiffness terms. 
The stiffness of the boom is modeled as 
fy, = A,u? + B, u?+C, u+ dy 

and the damping as 

fp, = A,U? + Bou + C 
For the present analysis, the damping is assumed to be linear 
and contains only the linear term corresponding to equivalent 
viscous damping. The nonlinear stiffness contribution due to 
the joints is modeled by the cubic term, A,u>. The linear term 


representing the boom stiffness distribution is determined 


from the experimental stiffness properties of the beam using 


the relation 





12ET 


Ke 2 


(74) 


The coefficient for the cubic term is treated as a parameter 
and as part of the parametric study, is varied as a function 


of the linear term as follows 
A, = C, + nC, (75) 


where n varies from t.1 to +.5. The equation of motion reduces 


to Duffing’s equation in the following form: 


mu + Bou + Aju? + Clu = Acoswt (76) 

The response of this nonlinear system to a forcing function is 
determined by approximating the derivatives as shown in the 
equation of motion. A solution based on step-by-step 
integration is used to generate the time response. A FORTRAN 
Program was developed based on the "linear acceleration 
method" and is presented in Appendix I. The recursion formulas 
used for the numerical integration is derived in the following 
analysis. 

in the linear acceleration method adopted here, the 
acceleration is approximated for a given step by the following 
relation: 


Ps Aud 
U(t) = Uy + art (77) 





Integration of equation (77) yields 
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= (78) 





and 


2 Ad, 23 
Uy, = Uy + Ut + us + ate (79) 


By using incremental quantities, Ap,, Au,, Al,, and Au, the 
computational algorithm 1s set up. Equation (79) may be solved 


for u, and equations (78) and (79) are be combined to give Al, 


as follows: 
= —6 - 6 ai 
BOP e pacts = sag he Shs (20) 
= 3 At 
a Maer 0a! Mala a 


Since equation (73) is satisfied at both t, and t,,,, it may be 


written as 
mA, + C,Au, + k,Au, = Ap, (82) 


Combining equations (80), (81), and (82) yields the 


incremental equation of motion 


KjAu, = Ap} (83) 
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Kl=K, -(z2) + 2 (84) 


i 
A. 





Ap: = Ap, + ig + 3C.)u, + [>= + cae] a, (85) 


Once Au, is determined from equation (83), Al, is obtainecé <ron 
equation (83) enc At, from equation (80). The updated values 


of u, vw, and t az step (+3) is computed Eroa 
=U, + Ad, (86) 


This step by step procedure described determines the 
response of the LACE spacecraft boom modeled as a nonlinear 
system. The numerical integration provides a time history of 
the response, namely deflection, u, velocity, v, and 
acceleration, a. 

Since experimental data for a 23-foot LACE-model boom is 
available from AEC-ABLE Engineering, the length of the boom 
analyzed is taken to be 23 feet. In the parametric study, the 
parameters varied are the excitation amplitude, frequency and 
she coefficient of the cubic stiffness term. The excitation 
amplitudes used gives deflections of 2.5%, 5%, 10%, 20%, and 


25% of the length of the spacecraft boom. The frequencies used 
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are 1/3, and 3w,n. These fxecuencies are chosen based on the 
exubic term of the stiffness that cenerates subharmonic and 
superharmonic response at one-third and three tizes the 
excitation frecuency, respectively. The coefficient of the 


cubic stiffness term is varied from 4.1 to +.5 of the linear 


The constant terms used in the equation of motion are as 
follows: 


= 3.02 lbs/iz 


a 


m= .1475 slugs 


a 
u 


-177 slugs/sec 
w, = 4.6 rad/sec 

where ¢=2/km{ - The damping factor{ was obtained from AEC- 
ABLE Engineering using the Logarithmic Decrement method. A 23 
foot boom mounted in its deployment cannister rigidly attached 
to a wall yields {7=.039 and a 23-foot boom rigidly mounted to 
a@ wall yields €=.001. The first case approximates the LACE 
spacecraft configuration. Both damping cases are considered in 


the present simulation. 


D. CHAOTIC VIBRATION ANALYSIS 

A comprehensive FORTRAN program called CHAOS was ieveloped 
by Sarigul-Klijn (Ref. 22] to analyze chaotic vibratiuns. The 
time history data from the numerical simulation program is 
used by the CHAOS program for the chaotic vibration analysis. 


The program CHAOS currently has 12 analysis options. They are 
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1. Plot the time history 

2. Plot onto toroidal phase space 

3. Take Poinceré section of toroida! phase space 
4- Plot onto 2-D phase space 

5. Plot onto 2-D Van Ger Pol plane 

6. Plot onto 3-D phase space 

7. ‘ake Poincaré section of 3-D phase space 

8. 3-D Poincaré section of a 4-D hyperspace 

9. Compute fourier power spectrum 

10. Compute statistics 

11. Compute Lyapunov exponent 

12. Compute fractal correlation d:mension 

As discussed earlier, the stiffness term is modeled as a 
combination of linear and cubic terms. The coefficient of the 
cubic term is varied as a percentage of the linear term from 
10% to 50%, and is seen to yield no significant changes in the 
response. Therefore, an average value of 30% is used in all 
subsequent numerical simulation. 

In order to focus the present study within the scope of 
this research, the study is divided into four specific groups. 
Each group is identified by excitation frequency and damping. 
In each group, the excitation amplitudes are 20.8 lbs, 41.7 
ibs, 83.4 lbs, 166.7 lbs, and 208.4 lbs respectively. This 
yields a deflection of 2.5%, 5%, 10%, 20%, and 25% of the 
length of the beam, respectively. Table VI contains the group 


identification. 
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Table VI GROUP IDENTIFICATION 


GROUP 





The initial analysis consists of studying four qualitative 
indicators of chaos. They are time history, fourier power 
spectrum, pseudo-phase plane and Poincare section 
respectively. Then, quantitative measures of the motion are 
computed. 

Figures 39-43 and 44-48 present the qualitative indicators 
for cases A and B respectively. The time history is initially 
transient and then becomes periodic in both cases. The fourier 
power spectrum indicates a broad-band frequency spectrum, with 
the number of spikes increasing with an increase in amplitude. 
The pseudo-phase plane plot indicates that the motion settles 
to a single limit cycle for case A and a double limit cycle 
for case B. The Poincaré section is very sparse and indicates 
a periodic system. The periodic time history also indicates a 
non~chaotic response. 

Figures 49-53 and 54-58 present the qualitative indicators 
of chaos for groups C and D. In all cases, the time history is 
seen to be non-periodic and the Fourier power spectrum is 


observed to be broad-band. In both cases, there appears to be 


two limit cycles in the pseudo-phase plot. The Poincaré 
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sections appear more dense than the previous groups and has 
the appearance of a strange attractor. These systems appear to 
be chaotic. 

To ascertain the quantitative nature of the response, the 
fractal correlation dimension is computed for each group and 
is shown in Table VII. Figures 59 and 60 show the fractal 
dimension for increasing phase space dimension for two cases 
of damping with F=20.8 and 0=1.53 respectively. The plots 
for the other cases are given in Appendix J. 

A non~integer fractal dimension is indicative of the 
presence of chaos in a system [Refs. 22,23,24,25]. This is 
consistent with the qualitative indicators of groups C and D. 
The inconsistency with groups A and B indicate that the 
fractal dimension by itself does not confirm the strangeness 
of the strange attractor as believed in the reported 
literature and must be used in conjunction with other 


indicators of chaos. Additional results pertaining to this 


parametric study are reported elsewhere (Ref. 26). 
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Figure 39. Qualitative Indicators Of Chaos. 
F=20.8 lbs, 0 =1.53 Hz, C=.177. 
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Figure 40. Qualitative Indicators Of Chaos. 
F=41.68 lbs, 0=1.53 Hz, C=.177 
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Figure 41. Qualitative Indicators Of Chaos. 
F=83.352 lbs, 0=1.53 Hz, C=.177 
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Figure 42. Qualitative Indicators Of Chaos. 
F=166.7 lbs, 0=1.53 Hz, C=.177 
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Figure 43. Qualitative Indicators Of Chaos. 
F=208.38 lbs, 0=1.53 Hz, C=.177 
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Figure 45. Qualitative Indicators Of Chaos. 
F=41.68 lbs, 0=13.8 Hz, C=.177 
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Figure 47. Qualitative Indicators Of Chaos. 
F=166.7 lbs, 0=13.8 Hz, C=.177 
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Figure 49. Qualitative Indicators Of Chaos. 
F=20.8 lbs, 1 =1.53 Hz, C=.0454 
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Figure 50. Qualitative Indicators Of Chaos. 
F=41.68 lbs, N=1.53 Hz, C=.0454 
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Figure 51. Qualitative Indicators Of Chaos. 
F=83.352 lbs, 0=1.53 Hz, C=.0454 
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Figure 52. Qualitative Indicators Of Chaos. 
F=166.7 lbs, 0=1.53 Hz, C=.0454 
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Figure 53. Qualitative Indicators Of Chaos. 
F=208.38 lbs, 0=1.53 Hz, C=.0454 
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Figure 59. Fractal Dimension 
F=20.8, N=1,.53, C=.177 
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Figure 60. Fractal Dimension. 


F=20.8, N=1.53, C=.0454 
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VIII. CONCLUSIONS AND SCOPE FOR FURTHER RESEARCH 


A. CONCLUSIONS 

A detailed finite-element model of the LACE spacecrazt has 
been developed using the finite element program GIFTS. A 
Gynanic analysis was performed on the model to determine the 
natural freauencies and mode shapes. This mode] will provide 
@ basis for comparison when actual orbital test data is 
avaiiable. 

Thermoelastic effects of the LACE spacecraft boom were 
investigated for a worst case temperature scenario. The 
results indicated that the deflections induced would have 
negligible effect on the systems dynamics and would not affect 
the on-orbit systems identification or induced thermal 
flutter. 

A preliminary study was conducted on a 23-foot section of 
the LACE spacecraft boom to investigate the possibility of 
chaotic vibrations occuring by varying the excitation 
amplitude and frequency. The system was modeled as a single- 
degree-of-freedom system with arbitrary nonlinear stiffness, 
and nonlinear damping excited by a sinusoidal function. A 
numerical integration program was developed to determine the 
system response using various excitation amplitudes and 


frequencies. Four qualitative methods, time history, fourier 
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power spectrum, pse.do-plane, and Poincaré section methods 
were used to evaluate chaotic vibrations. In addition, 
quantitative analysis included computation of correlation 
fractal dimension. It was evidenced that a combination of ail 
appoaches was needed to determine the nature of vibrations. 
Chaotic vibrations were present for the low-damped system 
while periodic response was indicated for highly damped 
systems. 

A non-integer correlation fractai dimension, which is 
believed to be a quantitative indication of chaos in a system, 
was evidenced for both periodic and non-periodic responses of 
the system. This important result, which is further being 
confirmed from analysis of flight test data (Ref. 21] shows 
that the fractal dimension by itself is not enough to 
characterize the attractors found in nonlinear dynamical 
systems. Modeling the “looseness" of joints for triangular 
space trusses by Moon & Li [Ref. 27}, who report interesting 
chaotic behavior, was brought to the attention of the author 
during the write-up of this report, which reveal the 
importance of understanding the nonlinear dynamic behavior of 


large flexible structures. 


B. SCOPE FOR FUTURE RESEARCH 
Future research may be beneficial in two main areas. 
The first is the development of a multi-body dynamics model 


of the LACE spacecraft. A structural model using GIFTS has 
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already been developed. The resuits generated may be used to 


Rapeereenen 


model the flexibility effects, while other effects introduced 
due to on-orbit, such as gravity gradients effects, ait 
atmospheric drag, etc., may be evaluated. These models may be oy , 
used as baselines for comparison when the orbital test data 
from LACE becomes available. This comparison will provide an 
opportunity for improvements and adequacy of these computer 
models. 

The availability of the orbital test data will provide 
information needed to conduct a more detailed study of chaotic 


vibrations on the fully deployed LACE spacecraft boom and 


s 


assess the sensitivity of the highly flexible boom structure 

to accidental or intentional excitation. Assessing the exact 

nature of the response of large flexible space structures is | 
important for missious that require accurate pointing 


accuracies and in suppression and control of vibraticns. 
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APPENDIX A 


$ ! COM FILS TO GENERATE SIMPLE MODEL OF LACE SPACECRAFT AND ‘ ‘ 
$ ! DERIVE SHE NATURAL FREQUENCIES. * & 
$ 


! FILENAME IS LACE.COM 
DEL LACE.*;* 


$ BEAMCS 

RACE 

OLB/LACECS 

QUrT : -. 
$ BOULKM 2 
LACE 

OLB/LACEBM 


ad 
> 
aQ 
is] 


$ GENERATES SOLID CIRCULAR BEAM USING PROCESSOR BEAMCS 
$ FILENAME IS LACECS.SRC 

CIRCS/1/.9585/ 

END 


$ SIMPLE FINITE ELEMENT MODEL OF THE LACE SPACECRAFT USING 
BULKM 
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$ FILENAME IS LACEBM.SRC Z 
$ ESTABLISH THE KEYPOINTS OF THE MODEL 


KPOINT/1/0,0,0/2/0,1800,0/3/0,~-900,0/ . 
KPOINT/4/0,0,1800/5/50,0,0/ a 


$ DEFINE MATERIAL PROPERTIES : 
ELMAT, 4/1/3.£3,8.E6,.29, .000014752/ 


; $ DEFINE BEAM CHARACTERISTICS 


S$ GENERATE THE BEAMS 
SLINE, 10/112/1,2,31/5/113/1,3,16/5/114/1,4,31/5/ 


END 


$ GENERATES POINT MASSES FOR THE SPACECRAFT 
$ FILENAME iS LACEBC.SRC 
MASS /MASSP/1/6.4078E/2/.0906/3/.0906/4/.517/ 


END F : 








APPENDIX B 


an i 5, 
bn nenthelecmtnmte: imme emnttepraareN, 


The following calculations model structural honeycomb 


oo. , 
panels as aluminum panels with equivalent thickness and xt ; 
f 
weight. The stiffness for a honeycomb panel is given as ; 
= _£th? 
2(1-v) 
and the equation for an aluminum plate is 
fe Bt 
12 (1 -v?) 
where 
E = 9.9 x 10° Ibs/in? 
v = .33 
h = core thickness of honmeycomb panel 
t = face thickness 
Throughout this appendix the following notation is used. 
p = density > t 


V = volume 


W = weight 


where subscripts h indicates honeycomb core and f refers to 


aluminum face skins. A panel size of 54" x 90" is assumed. 


A. PRIMARY STRUCTURE 
h= .9 
-05 


u 


t 































1. Determine Stiffness of Honeycomb Panel 


inn as = 
. D= 15 .9x40°) (.05) (.9)" _ 2.25x10° lbs./in 


2(4-.33)? 


2. Determine Equivalent Aluminum Plate Thicknessu 


3 6) £3 
p= — Fe 2.250 «105 = 12:9 420°) 


12(1-v*) 2(1- .33?) 
t = .624 in 


3. Determine Honeycomb Panel Mass 
p,* = 6 lbs/ft? = 3.47 x 10° lbs/in® 
p,* = 9.81 x 107% lbs/in? 
V, = +9 x 54 x 96 = 4665.6 in? 
Vet = «1 x 54 x 96 = 518.4 in? 
W, = 4665.6 x 3.47 x 10° = 16.2 lbs 
518.4 x 9.81 x 107% = 50.8 lbs 


= 
* 
u 


seasterik items remain constant throughout the 
calculations 


Total weight = 67 lbs \ 


67 Ibs 


Total mass => ————=——_ 
386.4 in/sec? 


= .1734 (Ibs-sec*)/in (GIFTS units) 


4. Detexmine Equivalent Density 


Volume of aluminum panel = .624x54x96 = 3234.8 in? 


Equivalent density = 2735 (Iba ~ sec’) /in = 5.364x10°5 (lbs ~sec*) /in* 
3234.8 in? 
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a SE RE COOL EF RE TE To Sow Soe 





B. SECONDARY STRUCTURE 
1. Fixed and Deployable Solar Array Substrate 
h= .4 
v 
t = .05 ‘ 


~ £9-9x10%) (.05) (.4)? - 4.4410! lbs-in 


D 
2 (1-~.337) 


a. Determine Equivalent Aluminum Thickness 


4.44x104 lbs-in = (9.9 £108) 6? 
12 (1-.337) 
t = .3634 in 


b. Determine Honeycomb Panel Mass 
3 


i 


V, = +4 % 54 x 96 = 2073.6 in 


W, = 2073.6 x 3.47 x 10% = 7.195 lbs 


n 


Total weight = 58 lbs 
otal mass = 58/386.4 = .15 ibs-sec*/in 
c. Determine Equivalent Density i 


Volume of aluminum panel = 1884 ind i 


Equivalent density = - aa = 7.96 x10°5 lbs-sec?/in* 





2. Deployable Sensor Panels 
h= .65 


t = .05 


2 
(9.9 20%) (.05) (65) - 4.173 x10° lbs-in 
2(1-.337) 


ps 














a. Determine Equivalent Aluminum Thickness 


(9.9 x10%) t> 


4.173 10° = - 
12 (1 - .33%) 


t = .502 in 


b. Determine Honeycomb Mass 


V, = -65 x 54 x 96 = 3369.6 in? 
W, = 3369.6 x 3.47 x 10? = 11.7 lbs 
Total weight = 62.5 lbs 


Total mass = .1617 lbs-sec*/in 


c. Determine Equivalent Density 


Volume of aluminum panel = .502 x 54 x 96 = 2602.4 in? 


Equivalent density = +2647). 6.22 x 10°5 (lbs-sec*) /in* 





2602.4 





APPENDIX C 


$ ! COM FILE TO GENERATE THE MAIN SPACECRAFT BODY 
$ ! FILENAME IS XLACE.COM 


DEL LACE.*;* 
$ BEAMCS 
LACE 
OLB/LACECS 
QUIT 
$ BULKM 
LACE 
OLB/LACEBM 
QUIT 
$ EDITM 
LACE 
OLB/LACEM 
QUIT 
$ LOADBC 
LACE 
MASS 

: OLB/MASS 
QUIT 





$ GENERATES THE BEAM CROSS SECTIONS FOR THE LACE SPACECRAFT 

$ GENERATES THE FILE USING THE PROCESSOR BEAMCS . \ 
$ SOURCE FILE NAME IS LACECS.SRC \ 
ANGLE, 11/1/2.875,3/.25, .25/0,1.938/90/ 

TBEAM, 10/2/2.75,3/.25, .25/0,1.938/ 

CHANNEL/3/7.75,3/.125,.125/ 

CIRCH,1/9/5.605,5.5/0,5.5525/ 


END 


$ GENERATES LACE SPACECRAFT MAIN STRUCTURE 
$ GENERATES THE FILE USING THE PROCESSOR BULKM : i 
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$ FILENAME IS LACEBM.SRC 


$ ELEMENT MATERIAL DEFINITIONS 

$ FRAME ELEMENTS 

ELMAT, 4/1/1.8E4,9.9E6,.33,2.5382E-4/ 

$ PRIMARY STRUCTURE 

ELMAT, 4/2/1.8E4,9.986,.33,5.3645-5/ 

$ SECONDARY STRUCTURE-FIXED AND DEPLOYABLE SOLAR ARRAY 
SUBSTRATE 

ELMAT, 4/3/1.8E4,9.956,.33,8.1915E-5/ 

$ SECONDARY STRUCTURE-DEPLOYABLE SENSOR PANELS 

ELMAT, 4/4/1.8E4,9.9E6,.33,6.22E-5/ 

$ SECONDARY STRUCTURE-DEPLOYABLE SENSOR ARMS AND BOTTOM PANEL 
ELMAT, 4/5/1.8E4,9.9E6,.33,2.5382E-4/ 

$ DEFINE ELEMENT THICKNESS 

$ PRIMARY STRUCTURE PANELS 

ETH,1/4/.624/ 

$ SECONDARY STRUCTURE PANELS- SOLAR ARRAY SUBSTRATE 

ETH, 1/5/.3634/ 

$ SECONDARY STRUCTURE PANELS-DEPLOYABLE SENSOR PANELS 

ETH, 1/6/.502/ 

$ SECONDARY STRUCTURE PANELS- DEPLOYABLE SENSOR ARMS AND 
BOTTOM 

$ PANEL 

ETH,1/7/.125/ 


$ ESTABLISH THE KEYPOINTS FOR THE PRIMARY STRUCTURE 


KPOINT/1/26,26,0/2/~26,26,0/3/26,-26,0/4/-26,-26,0/ 
KPOINT/5/26,26,41.75/6/-26,26,41.75/7/26,-26,41.75/8/-26,-26 
741.75/ 


ard 











KPOINT/9/26,26,68.375/10/-26,26,68.375/11/26,-26,68.375/ ae 
KPOINT/ 12/-26,-~26,68.375/13/26,26,96/ | 
KPOINT/14/-26,26,96/15/26,-26,96/16/-26,-26,96/ an 
$ ESTABLISH KEYPOINTS FOR THE CIRCULAR ARCS 

$ SIDE C+Y 4 
KPOINT/17/13,26,68.375/18/13,26,76.5/19/13,26,87.5/20/13,26, : 
96/ 

KPOINT/21/18.5,26,82/22/7.5,26,82/ 

$ SIDE C-Y 


KPOINT/23/-13,-26,68.375/24/-13,~-26,76.5/25/-13,-26,87.5/ ' 
KPOINT/26/-13,~-26,96/27/-7.5,-26,82/28/-18.5,-26,82/ 


$ SIDE +2 


KPOINT/29/26,0,96/30/5.5,0,96/31/-5.5,0,96/32/-26,0,96/ 
KPOINT/33/0,-5.5,96/34/0,5.5,96/ 


$ GENERATE VERTICAL ANGLE BEAM MEMBERS FOR C~SECTION PANELS 
LETY/BEAM2/1,1/ 

SLINE, 10/L913/9,13,6/10/L1014/10,14,6/9/L1216/12,16,6/10/ | 
SLINE,10/L1115/11,15,6/9/ 

$ GENERATE VERTICAL ANGLE BEAM MEMBERS FOR B-SECTION PANELS ‘ 


SLINE,10/L59/5,9,3/6/L610/6,10,3/5/L812/8,12,3/7/L711/7,11,3 
/5/ 


$ GENERATE VERTICAL ANGLE BEAM MEMBERS FOR A-SECTION PANELS 
SLINE,10/L15/1,5,3/2/L26/2,6,3/1/L48/4 ,8,3/3/L37/3,7,3/1/ 

$ GENERATE HORIZONTAL ANGLE BEAMS FOR TOP AND BOTTOM MEMBERS 

$ BOTTOM MEMBERS 

SLINE/L12/1,2,6/L34/3,4,6/L13/1,3,5/L24/2,4,5/ 

$ TOP MEMBERS 


SLINE/L1320/13,20,2/L1420/14,20,5/L1432/14, 32, 3/L1632/16,32,3/ 
SLINE/L1626/16,26,2/L1526/15,26,5/L1329/13, 29,3/L1529/15,29,3/ 


$ GENERATE HORIZONTAL T-SECTION BEAM MEMBERS 
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$ B-SECTION PANELS 


LETY/BEAM2/1,2/ 7 
SLINE, 10/L56/5, 6, 6/7/£78/7 ,8, 6/5/257/5,7,5/6/258/6,8,5/3/ 
$ C-SECTION PANEL 

SEINE, 10/L917/9, 17, 2/13/L1017/20,27,5/24/£2012/210,12,5/9/ 
SLINE, 10/L1223/12,23,2/10/L1123/21,23,5/9/£921/9 ,23,5/20/ 


$ GENERATE LINES ASSOCIATED WITH CIRCULAR ARCS 


$ GENERATE CIRCULAR ARCS 

$ C+¥ SIDE PANEL 

CARC/C1819/18,22,19,4/C1928/29,21,18,4/ 

$ C-¥ SIDE PANEL 

CARC/C1819/18,22,19,4/C1918/19,23,18,4/ 

$ C-¥ SIDE PANEL 

CARC/C2425/24,28,25,4/C2524/25,27,245,4/ 

$ TOP PANEL +2 

CARC/C3031/30,33,31,4/C3130/31,34,30,4/ 

$ GENERATE COMPLINES 
COMPLINE/L1314/L1320,L1420/L1416/L1432,L1632/1-1315/L1329,L15 
COMPLINE /L1516/11526,11626/%1720/11718,C1619,1.1920/ 
COMPLINE/L2017/L1920,C1918,1L1718/L2326/L2324 ,C2425,L2526/ 
COMPLINE/L2326/L2324 ,C2425,L2526/L2623/L2526,C2524,L2324/ 
COMPLINE/L910/L917,L1017/L1112/L1123,L1223/ 
COMPLINE/L2932/L2930,C3031,1L3132/L3229/L3132,C3130,L2930/ 
$ GENERATE GRID ELEMENTS 

GETY/QB4/2,4/ 

$ A-PANELS 


GRID4/A+Y/L12,L56,L15,L26/A-X/L24,L68,L26,L48/ 
GRID4/A-Y/L34,L73,1L48,L37/A+X/L13,L57,L15,L37/ 
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$ B-PENELS 


GRID4 /3B+F/L56,19120,259,3610/3-%/L68,41012,2610,L832/ 
GRIDS /3-¥/L78 ,L1212,2822,£721/3+X/L57,L911,L59, 2722 


$ GENERETE C-DANELS 

GRIDS /C+¥L/1917 ,£2017 ,51320,3913/C+ER/L1017 ,£1014, 52426,22720/ 
GRID4 /C-¥3/L1223 ,L2325,51626 ,L1216/C-2R/L1223,11115,41526,1L2 
Gath¢ /C4X/L913,11315,11115,1911/C-X/L1014,11¢16,11216,12022/ 
$ GENERATE TOP PANEL 

GRIDS /ZL/L1529 ,L1516 ,£1632,22932/ZR/L1329 ,L3229 ,L1432,L1314/ 
$ GENERATE RF PANEL 

GRiD4/RF/L56,1L68,578,L57/ 


$ GENERATE BOTTOM PANEL 


GETY/Q34/5,7/ 
GRID4 /BOTTOM/L12,124,L34,L13/ 


$ GENERATE KEYPOINTS FOR SENSOR PANELS AND SENSOR ARMS 


KPOINT/35/19,26,0/36/19,75.25,0/37/-19,26,0/38/~19,75.25,0/ 
KPOINT/39/-26,19,0/40/-75.25,19,0/41/-26,-19,0/42/-75.25,-19 
0/ 
KPOINT/43/-19,-26,0/44/-19,-75.25, 0/45/19, -26,0/46/19,~75.25 
10/ 
KPOINT/47/26,-19,0/48/75.25,~-19,0/49/26,19,0/50/75.25,19,0/ 
KPOINT/51/79.03,79.03,0/52/-79 .03,79.03,0/53/-79.03,-79.03,0/ 
KPOINT/54/79.03,-79.03,0/ 


$ GENERATE LINES FOR SENSOR PANELS 

SLINE/L3536/35, 36, 5/L3738/37 , 38 ,5/L3537/35, 37, 4/L3638/36,38,4/ 
SLINE/L4344/43,44,5/L4546/45 ,46,5/L4345/43,45,4/L4446/44,46,4/ 
SLINE/L4748/47,48,5/L4950/49 ,50,5/L4749/47,49,4/L4850/48,50,4/ 
$ GENERATE GRID ELEMENTS 

GETY/QB4/4,6/ 


GRID4/+¥/L3536 ,L3738,1L3537,L3638/-X/L3940,L4142,1L3941,14042/ 
GRID4/-Y/L4344,L4546,L4345,L4446/+X/L4748,L4950,L4749,L4850/ 


$ GENERATE CHANNEL BEAM ELEMENTS FOR SENSOR ARMS 








LETY/SEAM2/3,i/ 
SLINZ, 10/1151/1,53,6/52/L252/2,52,6/51/2453/4,53,6/54/ 
SLINE, 10/L354/3,54,6/53/ 


$ GENERATE KEXYPOINTS FOR SOLAR PANEL FOLDOUTS 


KPOINT/55/24,26,96/56/24,68.5,96/57/-24,26,96/58/-24 ,68.5,96/ 
KPOINT/59/-26,24,96/60/-68.5,24,96/6€1/-26,-24,96/62/-68.5,-2 
4,9C/ 
KPOINT/63/-24,-26,96/64/-24,-68.5,96/65/24,-26,96/66/24, -68. 
5,98/ 
KPOINT/67/26, -24,96/68/68.5,-24, 96/69/26,24,96/70/68.5,24,96/ 


$ GENERATE LINES FOR SOLAR PANEL FOLDOUTS 

SLINE/L5556/55,56, 3/L5758/57 ,58,3/L5557/55,57,3/L5658/56,58,3/ 
SLINE/L5960/59,60,3/L6162/61,62,3/L5961/59, 61,3/L6062/60,62,3/ 
SLINE/L6364/63,64,3/L6566/65,66,3/L6365/63,65,3/L6466/64,66,3/ 
SLINE/L6768/67 , 68, 3/L6970/69 ,70,3/L6769/67, 69, 3/L6870/68,70,3/ 


$ GENERATE GRID ELEMENTS 





GETY/QB4/3,5/ 
GRID4/D+¥/L5556 ,L5758,L5557 , L5658/D-X/L5960 ,L6162,L5961,L6062/ : 
GRID4/D-Y/L6364 , L6566, L6365, L6466 /D+X/L67 68 ,L6970,L6769 , L6870/ 
$ GENERATE KEYPOINTS FOR SOLAR FIXED PANELS | 
$ A+Y PANELS j | 


KPOINT/71/-26,26,20/72/-26,26,46/73/-14,26,46/74/-14,26,60/ 
KPOINT/75/14,26,60/76/14,26,46/77/26,26,46/78/26,26,20/ 


$ A-Y | 
KPOINT/79/26,-26,20/80/26,-26,46/81/14,-26,46/82/14,-26,60/ 

a a a a a i a a a , 

$ C+X 
KPOINT/87/26,26,30/88/26,26,45/89/26,14,45/90/26,14,56.01/ 
KPOINT/91/26,26,56.01/92/26,26,88.5/93/26,-26,88.5/94/26,-26 
SoOURT 95/3 654K, 8 6s01/96/20/48-48 97/26 206 48/90/26 38) S07 
$ C-X 

BE te alan can Baan engage ae ae 
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: KPOINT/103/-26,26,56.01/104/-26,26,88.5/105/-26,-26,88.5/ 
KPOINT/106/-26,-26,55.01/107/-26,-14,565.01/108/-26,-14,45/ 
KPOINT/109/-26,-26,45/110/-26,-26,30/ 


END * , 


$ EDIT THE MODEL-GENERATE VERTICAL ANGLE BEAM ELEMENTS 
$ GENERATES THE FILE USING THE PROCESSOR EDITM 

$ FILE NAME IS LACEM.SRC 

PTRM/1/PTRTH/1 

$ GENERATE TOP ANGLE MEMBERS 
BEAM2/13,55,29/20,55,29/20,57,32/14,57,32/ 
BEAM2/14,59,20/32,59,31/32,61,31/16,61,26/ 
BEAM2/16,63,32/26,63,32/26,65,32/15,65,29/ 
BEAM2/15,67,26/29,67,26/29,69,30/13,69,20/ 

$ GENERATE BOTTOM ANGLE MEMBERS 
BEAM2/1,35,3/35,37,3/2,37,4/ 
BEAM2/2,39,1/39,41,1/4,41,3/ 
BEAM2/4,43,2/43,45,2/3,45,1/ 
BEAM2/3,47,4/47,49,4/1,49,2/ 
BEAM2/419,422,1/411,414,1/135,139,1/138,142,1/ . ‘ 


$ GENERATE HOLLOW CIRCULAR BEAM CROSS SECTION FOR CYLINDRICAL 
$ CANNISTERS 


PTRTH/9/ 
BEAM2/19, 364,9/25,352,9/31,405,9/ 
END 


GENERATES POINT AND GRID MASSES FOR THE LACE SPACECRAFT 
GENERATES THE FILE USING THE PROCESSOR LOADBC 
SOURCE FILE NAME IS MASS.SRC 


GENERATE MASS DISTRIBUTIONS FOR A-PANELS 


nn wm Ww mn 


MASS LOADING FOR A-Y PANEL 





MASSP/8/6.31E~-3/166/1.54E-2/165/1.65E-2/ 
MASSP/164/1.65E-2/163/1.19E-2/7/2.87E-3/ 
MASSP/133/9.76E-3/325/2.48E-2/324/2.25E-2/ 
MASSP/323/2.25E-2/322/1.652-2/134/7.4E-3/ 
MASSP/4/6.31E~3/142/1.23E-2/141/8.88E-3/ 
MASSP/140/8.885-3/139/7.4E-3/3/7.4E-3/ 


$ MASS LOADING FOR A+Y PANEL 


MASSP/5/2.87E-3/159/1.37E-2/160/2.03E-2/ 
MASSP/161/1.61E-2/162/9.47£-3/6/2.87E-3/ 
MASSP/131/5.27E-3/315/1.61E-2/316/2.69E-2/ 
MASSP/317/2.93E-2/318/1.61E-2/132/2.87E-3/ 
MASSP/1/5.27E-3/135/5 .27E-3/136/9 .47E-3/ 
MASSP/137/1.61E-2/138/9.47E-3/2/2.87E-3/ 


$ MASS LOADING FOR A-X PANEL 


MASSP/6/7.01E-3/170/1.86E-2/171/3.07E-2/ 
MASSP/172/2.36E-2/8/7.97E-3/ 
MASSP/132/2.55E-2/319/5.57E-2/320/7.22E-2/ 
MASSP/321/4.66E-2/133/7.975-3/ 
MASSP/2/2.2E-2/146/4.05E-2/147/4.48E~-2/ 
MASSP/148/2.64E-2/4/3.44E~3/ 


$ MASS LOADING FOR A+X PANEL 
MASSP/7/3.44E-3/169/1.88E-2/168/2.3E-2/ 
MASSP/167/7.66E-3/5/3.44E-3/ 
MASSP/134/1.32E-2/328/4.7E-2/ 
MASSP/327/6.0E-2/326/3.89E-2/131/1.62E-2/ 
MASSP/3/1.32E-2/145/3.16E-2/144/4.03E-2/ 
MASSP/143/3.46E-2/1/1.62E-2/ 

$ GENERATES MASS DISTRIBUTIONS FOR B-PANELS 
$ MASS LOADING FOR B+Y PANEL 
MASSP/5/8.4E~3/127/8.4E-3/329/8.4E-3/159/8.4E-3/ 
MASSP/160/1.92E-3/330/1.92E-3/161/3.4E-2/ 
MASSP/331/3.4E-2/162/1.59E-2/332/1.59E-2/ 

$ MASS LOADING FOR B-X PANEL 
MASSP/170/7.47E=3/171/7 .47E-3/172/7.47E-3/ 
MASSP/333/7.47E~3/335/7.47E-3/ 
MASSP/334/8.9E-3/177/2.16E-3/ 


$ MASS LOADING FOR B+X PANEL 
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MASSP/11/4.27E~3/184/8.83E-3/183/9.81E-3/ 
MASSP/182/1.7E-2/9/1.24E-2/130/4.27E-3/ 
MASSP/342/1.34E-2/341/1.89E-2/340/3.48E-2/ 
MASSP/127/2.57E-2/7/6.91E-4/169/5.23E-3/ 
MASSP/168/9.81E-3/167/1.79E-2/5/1.39E-2/ 


$ GENERATES MASS DISTRIBUTIONS FOR C-PANELS 
$ MASS LOADING FOR C+Y PANEL 


MASSP/344/3.96E-3/345/3 . 96E-3/347/3 .96E-3/ 
MASSP/348/3 .96E-3/174/3.96E-3/175/3.96E~3/ 


$ MASS LOADING FOR C-Y PANEL 


MASSP/181/4.84E-3/357/4.84E-3/360/9 .68E-3/ 
MASSP/363/4.84E-3/366/4.84E-3/ 
MASSP/180/4 .84E-3/356/4 .84E-3/359/9.68E-3/ 
MASSP/362/4.84E-3/365/4.84E-3/ 


$ MASS LOADING FOR C-X PANEL 


MASSP/379/8 .67E-3/380/8 .67E~3/383/8.67E-3/ 
MASSP/384/8 .67E~3/387/8 .67E-3/388/8.67E-3/ 


$ MASS LOADING FOR C+X PANEL 


MASSP/158/3.35E-3/378/3.35E-3/368/3.08E-3/ 
MASSP/369/3 .08E-3/372/3.08E-3/373/3.08E-3/ 
MASSP/371/3.81E-3/ 


$ GENERATES MASS DISTRIBUTIONS FOR THE TOP 


MASSP/16/5 .29E-3/153/4.21E-3/32/3.56E-3/ 
MASSP/152/9.7E-4/14/2.06E-3/15/2.06E-3/ 
MASSP/158/9.7E-4/29/9.7E-4/157/9.7E-4/ 
MASSP/13/2.06E-3/26/5.71E-3/156/3.77E-3/ 
MASSP/155/3.77E-3/154/2.47E-3/149/2.47E-3/ 
MASSP/150/3.44E-3/151/3.44E-3/20/2.47E-3/ 
MASSP/394/4.21E-3/393/2.26E-3/392/2.26E-3/ 
MASSP/391/9.7E-4/30/9.7E-4/193/9.7E-4/ 
MASSP/194/9.7E-4/31/9.7E-4/195/9.7E-4/ 
MASSP/196/9.7E-4/398/9.7E-4/397/1.94E-3/ 
MASSP/396/1.94E-3/395/9.7E-4/ 


PANEL 


$ GENERATE MASS DISTRIBUTIONS FOR THE RF PANEL 


MASSP/8/2.16E-3/172/2.16E-3/171/2.52E-3/ 
MASSP/170/2.16E-3/6/2.16E-3/166/5.49E-3/ 
MASSP/410/1.29E-2/406/8.34E-3/402/2.16E-3/ 
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MASSP/162/2.16E-3/165/6.0E-3/409/1.18E-2/ 
MASSP/405/1.66E-2/401/1.07E-2/161/1.07E-2/ 
MASSP/164/6.0E-3/408/9 .88E-3/404/2.32E-2/ 

MASSP/400/1.93E-2/160/1.93E-2/163/5.49E-3/ 
MASSP/407/9.37E-3/403/1.46E-2/399/1.07E-2/ 
MASSP/159/1.07E-2/7/2.56E-3/169/6.03E-3/ ‘ J 
MASSP/168/6.32E-3/167/2.16E-3/5/2.16E-3/ 


$ GENERATES MASS DISTRIBUTION ON THE BOTTOM PANEL 


MASSP/4/5.85E-3/148/9.65E-3/147/1.06E-2/ 
MASSP/146/8.36E-3/2/5.85E-3/142/3.19E-2/ 
MASSP/422/3.77E-2/418/1.73E-2/414/2.23E-2/ 
MASSP/138/1.28E-2/141/2.83E-2/421/5.5E~2/ 
MASSP/417/6.07E-2/413/4.34E-2/137/1.15E-2/ 
MASSP/140/2.57E-2/420/5.8E-2/416/6.21E~2/ 
MASSP/412/3.39E-2/136/5.91E-3/139/2.94E-2/ 
MASSP/419/3.78E-2/415/1.9E-2/411/1.56E~2/ 
MASSP/135/7.14E-3/3/5.85E-3/145/7.14E-3/ 
MASSP/144/8.62E-3/143/1.1E~2/1/5.85E-3/ 


$ GENERATES MASS DISTRIBUTIONS ON THE SENSOR PANELS AND SENSOR 
ARMS 


$ GENERATE MASS LOADINGS FOR THE SENSOR ARMS 


MASSP/1/8.16E-3/283/7 .89E-3/284/1.12E-2/ 

MASSP/285/7 .89E-3/286/8.16E-3/51/7 .89E-3/ 

MASSP/2/8.16E-3/287/7 .89E-3/288/1.12E-2/ ‘ ‘ 
MASSP/289/7 .89E-3/290/8.16E-3/52/7 .89E~3/ 
MASSP/3/8.16E-3/295/7.89E-3/296/1.12E-2/ 
MASSP/297/7.89E-3/298/8.16E-3/54/7.89E-3/ 
MASSP/4/8.16E-3/291/7.89E-3/292/1.12E~2/ \ 
MASSP/293/7.89E-3/294/8.16E-3/53/7.89E-3/ 


$ GENERATE MASS LOADINGS FOR THE SENSOR PANELS 
$ GENERATE MASS LOADINGS FOR +Y PANEL 


MASSP/37/7.83E~-3/246/6.47E-3/247/5 .56E-3/ 
MASSP/248/5.63E-3/38/6.02E-3/ 
MASSP/250/4.66E-3/426/1.22E-2/ 
MASSP/427/1.22E-2/428/7.18E-3/252/4.66E-3/ 
MASSP/249/4.66E-3/423/1.22E~2/424/1.22E-2/ 
MASSP/425/7.18E~3/251/4.66E-3/ 
MASSP/35/6.02E-3/243/6.47E-3/ 
MASSP/244/5.56E~3/245/5 .63E-3/36/6.02E-3/ 


$ GENERATE MASS LOADINGS FOR -X PANEL 











































MASSP/41/3.49E-3/256/4 .66E-3/257/3.49E-3/ 
MASSP/258/4.14E-3/42/5.432-3/ 
MASSP/260/1.87E-2/431/1.93E-2/ 
MASSP/432/6.66E-3/433/6.28E-3/262/4.08E-3/ 
MASSP/259/1.87E-2/429/2.11E-2/430/6.66E-3/ 
MASSP/431/6.28E-3/261/4 . 08E-3/39/5.3E-3/253/2.85E-3/ 
MASSP/253/2.85E-3/254/3.49E-3/255/4.14E-3/40/5.43E-3/ 


$ GENERATE MASS LOADINGS FOR -Y PANEL i | 


MASSP/44/6.3E-3/265/5.91E-~3/264/5.85£-3/ 
MASSP/263/6.75E-3/43/8.11E-3/271/4.94E-3/ 
MASSP/437/7.47E-3/436/1.25E-2/435/1.25E-2/ 
MASSP/269/4.94E-3/272/4.94E-3/440/7 .47E-3/ 
MASSP/439/1.25E-2/438/1.25E-2/270/4.94E-3/ 
MASSP/46/6.3E-3/268/5.91E~3/267/5.85E-3/ 
MASSP/266/8.57E-3/45/6.3E~3/ 


$ GENERALE MASS LOADINGS FOR +X PANEL 


MASSP/48/5.14E-3/275/4.79E-3/274/3.45E-3/ 
MASSP/273/2.55E-3/47/4.79E-3/281/3.78E-3/ 
MASSP/443/5.98E-3/442/8.18E-3/441/3.78E-3/ 
MASSP/279/4.79E-3/282/3.78E-3/446/5.07E-3/ 
MASSP/445/8.18E-3/444/3.71E-3/280/3.78E-3/ 
MASSP/50/5.14E-3/278/4.79E-3/277/3.84E-3/ 
MASSP/276/2.48E-3/49/3.78E-3/ 


$ GENERATES GRID MASS FOR SOLAR PANELS 
$ SURFACE MOUNT SOLAR PANELS 


MASSG/A+Y/1.531815E-5,1.531815E-5,1.531815E~5,1.531815E-5/ 
MASSG/A-Y/1.5318155-5,1.£31815E-5,1.531815E-5,1.531815E-S/ 
MASSG/B+Y/2.402E-5,2.402E-5,2.402E-5,2.402E-5/ 
MASSG/B-Y/2.402E-5,2.402E-5,2.402E-5,2.402E~-5/ 
MASSG/A+X/1.196E-5,1.196E-5,1.196E-5,1.196E-S/ 
MASfS3/A-X/1.196E-5,1.196E-5,1.196E-5,1.196E-5/ 
MASSG/B+X/1.8755E-5,1.8755E-5,1.8755E-5,1.8755E-5/ 
MASSG/B-X/1.8755E-5,1.8755E-5,1.8755E-5,1.8755E-5/ 
MASSG/C+X/1.8076E-5,1.8076E-5,1.8076E-5,1.8076E-5/ 
MASSG/C-X/1.8076E-5,1.8676E-5,1.8076E-5,1.8076E-5/ 


Se Ne eee eNom MeReetAmannte A wh waren a He 
ea a Eee 


$ DEPLOYED SOLAR PANELS 


MASSG/D+X/6.343E-6,6.343E-6,6.343E-6,6.343E-6/ 
MASSG/D-X/6.343E-6,6.343E-6,6.343E-6,6.343E-6/ 
MASSG/D+Y/6.343E-6,6.343E-6,6.343E-6,6.343E-6/ 
MASSG/D-Y/6.343E-6,6.343E-6,6.343E-6,6.343E-6/ 


$ DETERMINES MISCELLANEOUS UNACCOUNTED MASS 
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$ DISTRIBUTED OVER SPACECRAFT PANELS 


$ A-PANELS - 
MASSG/A+K/1.46915E-5,1.4691E-5,1.4691E-5,1.4691E~5/ 
MASSG/A-X/1.4691E-5,1.4691E-5,1.46915-5,1.4691E-5/ ' ‘ 
MASSG/A+Y/1.4691E-5,1.4691E-5,1.4691E-5,1.4691E-5/ : . 
MASSG/A+Y/1.4691E-5,1.4691E-5,1.46915-5,1.4691E~5/ : 
$ B-PANELS 
MASSG/B+X/2.303662E-5,2.303662E-5,2.303662E-5,2.303662E-5/ 
MASSG/8-X/2.303662E-5,2.303662E-5,2.303662E-5,2.303662E-5/ 
MASSG/B+Y/2.303662E-5,2.303662E-5,2.303662E-8,2.3036625-5/ 
MASSG/B-X/2.303662E-5,2.303662E-5,2.303662E-5,2.303662E-5/ 
$ C-PANELS 


MASSG/C+X/2.2202728-5,2.220272E-5,2.220272E-5,2.220272E-5/ 
MASSG/C-X%/2.220272E-5,2.220272E-5,2.220272E-5,2.220272E-5/ 


$ TOP PANEL 


MASSG/ZL/5.897597E-6,5.897597E-6,5.897597E-6,5.897597E-6/ 
MASSG/2R/5.897597E-6,5.897597E-6,5.897597E-6,5.897597E-6/ 


$ RF PANEL 
MASSG/RF/1.226228E-5,1.226228E-5,1.226228E-5,1.226228E-5/ 
$ BOTTOM PANEL 


MASSG/BOTTOM/1.179519E-5,1.179519E-5,1.179519E-5,1.179519E-5/ 


END { 
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APPENDIZ E 


5 
ODE OF THE SPACECRAFT 





OL3/TRUSCS 

Quit? 

$ BULXM 

TRUS150 

OL3/TRS1503u 

QUIT 

$ LORDBC 

ore 
OL3/TRUSBC 

OLB/TRSMASI5 


TRUSi50 

6 
1,2,3,4,5,6 
OL3/TRUSDEF 
QUIT 

$ ADSTIF 
TRUS150 

$ STASS 
TRUSI5O 

$ DECOM 
TRUS150 

$ REDCS 
TRUS150 


! FILENAME IS TRUSCS.SRC 


0’ TRUSS RS & SUBSTRUCTURE TO BE 


! GENERATES CROSS SECTION TYPE AND DIMENSION OF BEAM ELEMENTS 


CIRCS/8/.075/END 
! BULKM FILE THAT GENERATES 150’ TRUSS 


! FILENAME IS TRUS150BM.SRC 


156 





ELERT, 4/6/6.0£4,8.056,.29,1.940994E-4/ 
KPOINT/1/0,0,4.763/2/0,8.25,0/3/0,0,-4.763/4/1800,0,4.763/ 
KPOE NT/5/1800, 8.25,0/6/1800,0,-4.763/7/90,0,50/ 

LETY /BEAM2/6,8/ 
SLINE,10/L14/1,4,301/7/L36/3,6,301/7/L25/2,5,301/8/L12/1,2,3 
/7/ 

SLINE, 10/L23/2,3,3/7/L13/1,3,3/7/%45/4,5,3/7/L56/5,6,3/7/ 
SLINE, 10/L46/4,6,3/7/ 

END 


! EDITM FILE TO ADD BATTENS AND DIAGONALS TO THE 150’ TRUSS 


ENAME IS TRSI50EM.SRC 


ct 


PTRTH/8/PTRIM/6/ 

$ VERTICAL MEMBERS 
BEAN2/9,307,7,349/157,455,7/158,307,7,149/306,455,7/ 

$ BOTTOM MEMBERS 

BEAM2/9,158,7,149/157,306,7/ 

$ DIAGONAL MEMBERS 

ETH,1/9/.001963/ 

PIRTH/9/ 
ROD2/9,159,148/156,306/10,158, 148/157, 305/9,308,148/156,455/ 
ROD2/10, 307, 148/157, 454/158, 308, 148/305, 455/159, 307,148/306, 
454/ 
ROD2/1,307/2,9/3,9/1,158/3,307/2,158/4,455/5,157/5,306/6,455/ 
ROD2/6,157/4,306/ 

END 

! GENERATES POINT MASSES FOR THE 150’ TRUSS 


! FILENAME IS TRSMAS15.SRC 


MASSP/1/2.71E-5/2/2.71E-5/3/2.71E-5/ 
MASSP/4/2.71E~-5/5/2.71E-5/6/2.71E-5/ 
MASSP/9,307,299/2.71E-5,2.71E-5/308, 606,299/2.71E~5,2.71E-5/ 
MASSP/6C7,905,299/2.71E-5,2.71E-5,/ 

END 





ee aE: ee ecm Hee 





POINT 
440 
441 
442 
443 
444 
445 
446 
447 
448 
449 
450 


mm OD eA ee OO ee OO ee OO 


MASS X 


- 234E-02 
. 650E-03 
. 305E-02 
- 985E-92 
. SBOE-03 
. 305E-02 
. S40E-03 
. 842E-92 
. 842E-02 
. 9B83E-03 
. 842E-02 


ee el oe eo 


MASS ¥ 


. 234£-02 
. §50E-03 
. 305E-92 
. O85E-02 
. 9B0E-03 
. 307E+00 
. 994E-02 
. 842E-02 
. 842E-02 
. 842E-02 
. 842E-02 





meee Oe Oe ee Oe 








MASS 2 


. 234E-02 
. §50E-03 
. 309E-92 
- 085E-02 
. S80E-93 
. 305E-02 
. §40E-03 
. 5B3E-03 
. 842E-02 
. 842E-02 
.§42E-02 


APPENDIX F 


$ ! COM FILE TO GENERATE AND CONNECT SPACECRAFT 


$ ! FILENAME IS TRUSS.COM 


OLB/TRUSCS 
Quit 

SBULKM 

LACE 
OLB/TRUSSBM 
Quit 

$ EDITM 
LACE 
OLB/TRUSSEM 
Quit 

$ LOADBC 
LACE 
OLB/TRUSMASS 
QUIT 


TRUSSES 


$ GENERATES THE FILE FOR THE THREE TRUSSES OF THE LACE 


$ SPACECRAFT USING THE PROCESSOR BULKM 


$ BULKM FILE TO GENERATE 150 FOOT TRUSS IN +Y DIRECTION 


$ ELEMENT MATERIAL 
ELMAT,4/6/6.0E4,8.026,.29,1.94099E~4/ 


KPOINT/451/8.237,26,79.25/452/8.237,1826,79.25/ 
KPOINT/453/13,26,87.5/454/13,1826,87.5/ 


KPOINT/455/17.763,26,79.25/456/17.763,1826,79.25/ 


KPOINT/457/67.763,926,79.25/ 
LETY/BEAM2/6,8/ 


SLINE, 10/LTR12/451, 452,301/457/LTR34/453,454,301/457/ 
SLINE, 10/LTR56/455, 456, 301/457/LTR13/451,453,3/457/ 
SLINE, 10/LTR35/453,455,3/457/LTR15/451,455,3/457/ 
SLINE, 10/LTR24/452,454,3/457/LTR46/454,456,3/457/ 


SLINE, 10/LTR26/452,456,3/457/ 


$ FILE TO GENERATE 150 FOOT TRUSS IN +2 DIRECTION 


159 


re 








’ 
i 
{ 
H 


$ ELEMENT MATERTAL 
ELMAT,4/6/6.054,8.056, .29,1.94099E-4/ 


KPOINT/458/2.75,4.763,96/460/-5.5,0,96/ 
KPOINT/462/2.75,-4.763,96/459/2.75,4.763,1896/ 
KPOINT/461/-5.5,0,1896/463/2.75,-4.763,1896/ : ’ 
KPOINT/470/52.75,0,996/ 


LETY/BEAN2/6,8/ 


SLINE, 10/LTR89 /458, 459 , 301/470/LTR6263/462,463,301/470/ 
SLINE,10/LTRO1/460, 461, 301/470/LTR5860/458, 460,3/470/ 
SLINE, 10/LTR6062/460, 462, 3/470/LTR5862/458,462,3/470/ 
SLINE, 10/LTR5961/459 , 461,3/470/LTR6163/461, 463,3/470/ 
SLINE, 10/LTR5963/459 , 463,3/470/ 


$ GENERATES A FILE TO CONNECT 75 FOOT -Y-TRUSS BOOMS TO 
$ MAIN SPACECRAFT 


KPOINT/464/-17.763,-26,79.25/ 
KPOINT/466/-13,-26,87.5/ 
KPOINT/468/-8.237,-26,73.25/ 
KPOINT/465/~-17.763,-926,79.25/ 
KPOINT/467/-13,-926,87.5/ 
KPOINT/469/-8.237,-926,79.25/ 
KPOINT/471/67.763,-463,79.25/ 


LETY/BEAM2/6,8/ 

SLINE, 10/L46869/468, 469, 151/471/246465/464,465,151/471/ 

SLINE, 10/L46667/466, 467, 151/471/L46567/465,467,3/471/ 

SLINE, 10/L46769/467,469,3/471/L46569/465,469,3/471/ 

SLINE, 10/L46466/464,466,3/471/146668/466,468,3/471/ : 
SLINE, 10/L46468/464,468,3/471/ 

END 

$ GENERATES BATTENS AND DIAGONALS FOR SPACECRAFT TRUSS 

$ USING THE PROCESSOR EDITM 

$ FILENAME IS TRUSSEM.SRC 

$ THIS FILE MERGES THE TRUSSES TO THE SPACECRAFT ' 
MERGEP/546/451/19/453/549/455/ | 


MERGEP/557/548/31/461/553/462/ 
MERGEP/550/464/25/466/553/468/ | 


160 





VERTICAL/BOTTOM MEMBERS 
EAN2/974,1273,457,299/1272,1571,457/ 
EAM2 /675,974,457,299/973,1272,457/ 


GENERATES BATTENS AND DIAGONALS FOR SPACECRAFT TRUSS 
USING THE PROCESSOR EDITM 


$ 
B 
B 
BEAM2/675,1273,457,299/973,1571,457 
$ 
$ 
$ 


FILENAME IS TRUSSEM.SRC 
$ THIS FILE MERGES THE TRUSSES TO THE SPACECRAFT 
MERGEP/546/451/19/453/549/455/ 
MERGEP/557/548/31/461/553/462/ 
MERGEP/550/464/25/466/553/468/ 


$ GENERATES VERTICAL, BOTTOM, DIAGONAL MEMBERS FOR 150 FT TRUSS 
$ IN THE +¥ DIRECTION 


PTRTH/8 
PTRM/6 


$ VERTICAL/BOTTOM MEMBERS 
BEAM2/974,1273,457,299/1272,1571,457/ 
BEAM2/675,974,457,299/973,1272,457/ 
BEAM2/675,1273,457,299/973,1571,457 

$ DIAGONAL MEMBERS 

ETH, 1/10/.001963/ 

PTRTH/10/ 

ROD2/455 ,675/451,1273/455,974/453,1273/453,675/451,974/ 
ROD2/975,1273,298/1272,1570/974,1274,298/1271,1571/ 
ROD2/675,975,298/972,1272/676,974,298/973,1271/ 
ROD2/675,1274,298/972,1571/676,1273,298/973,1570/ 
ROD2/456 ,973/452,1571/456,1272/454,1571/452,1272/454,973/ 


$ GENERATES VERTICAL, BOTTOM, DIAGONAL MEMBERS FOR 150 FT TRUSS 
$ IN THE +2 DIRECTION 


PTRTH/8 
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Ye OOAL Se eh ae ennai nena 


4 














PTRM/6 

$ VERTICAL/BOTTOM MEMBERS 
BEAM2/1578, 1877, 464,299/1876,2175,464/ 
BEAM2/1578,2176,464,299/1876,2474,464/ 
BEAM2/1877,2176,464,299/2175,2474,464 

$ DIAGONAL MEMBERS 

ETH,1/10/.001963/ 

PTRTH/10/ 
ROD2/463,1876/459,2175/461,1876/459,2474/461,2175/463, 2474/ 
ROD2/1877, 2177, 298/2174,2474/1878, 2176, 298/2175,2473/ 
ROD2/1579,1877,298/1876,2174/1578,1878,298/1875,2175/ 
ROD2/1579,2176,298/1876,2473/1578,2177, 298/1875,2474/ 
ROD2/462,2176/462,1578/460,1877/460,1578/458,2176/458,1877/ 


$ GENERATES VERTICAL, BOTTOM, DIAGONAL MEMBERS FOR 75 FT TRUSS 
$ IN THE -Y DIRECTION 


PTRTH/8 
PTRM/6 


$ VERTICAL/BOTTOM MEMBERS 

BEAM2 /2630,2779,470,149/2778,2927,470/ 
BEAM2/2481,2779,470,149/2629,2927,470/ 

BEAM2 /2481, 2630,470,149/2629,2778,470/ 

$ DIAGONAL MEMBERS 

ETH,1/10/.001963/ 

PTRTH/10/ 

ROD2/464,2779/466,2630/464,2481/468, 2630/468,2779/466,2481/ 
ROD2/2630,2780,148/2777 , 2927/2631,2779,148/2778, 2926/ 
ROD2/2482,2630,148/2629, 2777/2481, 2631, 148/2628,2778/ 
ROD2/2481,2780,148/2628, 2927/2482,2779,148/2629,2926/ 
ROD2/465,2927/467,2778/465,2629/469,2778/469, 2927/467, 2629/ 


END 


$ GENERATES MASS FOR THE COMPLETE TRUSS MODEL 


$ GENERATES MASS FOR THE 150 FOOT TRUSS IN +Y DIRECTION 























MASSP/451/8.8846E-5/452/8.8846E-5/453/8 .8846E-5/ 
MASSP/454/8.8846E-5/455/8.8846E-5/456/8.8846E-5/ 
MASSP/1273,1571,299/8.8846E-5,8.8846E-5/ 
MASSP/675,973,299/8.8846E-5,8.8846E-5/ 
MASSP/974,1272,299/8.8846E-5,8.8846E-5/ 


$ GENERATES MASS FOR THE 150 FOOT TRUSS IN +Z DIRECTION 


MASSP/464/8.8846E-5/465/8.8846E-5/466/8.8846E-5/ i 
MASSP/467/8.8846E-5/468/8 .8846E-5/469/8.8846E-5/ ‘ 
MASSP/2779,2927,149/8.8846E-5,8.8846E-5/ 
MASSP/2630,2778,149/8.8846E-5,8.8846E-5/ 
MASSP/2481,2629,149/8.8846E-5,8.8846E-5/ 


$ GENERATES MASS FOR THE 75 FOOT TRUSS IN -Y DIRECTION 


MASSP/458/8.8846E-5/459/8.8846E-5/460/8.8846E-5/ 
MASSP/461/8.8846E-5/462/8.8846E-5/463/8.8846E-5/ 
MASSP/1578,1876,299/8.8846E-5,8.8846E-5/ 
MASSP/1877,2175,299/8.8846E-5,8.8846E-5/ 
MASSP/2176,2474,299/8.8846E-5,8.8846E-5/ 


$ GENERATES TIPMASS FOR EACH TRUSS 
MASSP/459/1.72533E-1/461/1.72533E-1/463/1.72533E-1/ 
MASSP/452/3.0193E-2/454/3.0193E-2/456/3.0193E-2/ 
MASSP/465/3.0193E-2/467/3.0193E-2/469/3.0193E-2/ 


END 





APPENDIX G 


Before the LACE spacecraft was in its fully deployed state, 
it underwent three different configurations. Figures G1 to G4 
show the natural frequencies and mode shapes of the LACE 
spacecraft in its initial configuration with the gravity 
gradient extended to 75 feet. Figures G5 to G8 show LACE with 
the gravity gradient extended to 150 feet. Figures G9 to G12 
show the LACE spacecraft in its final configuration before 


full deployment with the lead and trail booms extended to 





119.5 feet. 
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Figure Gl. Mode i, u 
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Figure G2. Mode 2, v 
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Figure G3. Mode 3, u 





Figure G4. Mode 4, uv = .04579 Hz, GIFTS COMPLEX MODEL 
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.01217 Hz, GIFTS COMPLEX MODEL 


Mode 1, u 


Figure G5. 
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Figure G6. Mode 2, vu = .01226 Hz, GIFTS COMPLEX MODEL 
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Figure G7. Mode 3, v = .03535 Hz, GIFTS COMPLEX MODEL 
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Figure G8. Mode 4, v = .03552 Hz, GIFTS COMPLEX MODEL 


Figure G9. Mode 1, u = .02322 Hz, GIFTS COMPLEX MODEL 
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Figure G10. Mode 2, vu = .04294 Hz, GIFTS COMPLEX MODEL 
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Figure G12. Mode 4, v = .1349 Hz, GIFTS COMPLEX MODEL 








APPENDIX H 


‘ : A. MAXIMUM TEMPERATURE FOR DIAGONALS 


442 Btu/sq ft hr 
e=.8 
as = 1.0 


e=.8 70 Btu/sq ft hr 
8 160 Btu/sq ft hr 


Figure 11. LACE Boom Perpendicular to the Sun 


The heat balance equation is 


a,Sa = eoT'D 


where s is solar flux intensity, a is surface area ando is 


Stefan Boltzman constant. 


Heat in = 442 x 1.0 x d + (160 x 1.0 + 70 x .8) x nd/2 





—— 

















781.3 d 
Heat out = 0.8 xndxo x T! i 
781.3 = 0.9 xn xo x T = 4.305 x 10° 


T = 652 R= 193 F 


B. MINIMUM TEMPERATURE FOR THE DIAGONAL 


Direct solar 
442 cos 90°= 0 Btu/sq ft hr 


Earth emitted 
70 Btu/sq ft hr 


Earth reflected 
160 Btu/sq ft hr 


vA 


Figure H2. LACE Beam Parallel to the Sun 


Fa = +256 


: ‘ F,. 
Reciprocity Fy),dA, = F,.4,A, 4 at S 








ee 2 F;.. 
Divide by da, oe Asa, tt Fa Eb 
Orr, 256, 
. e- . 
Substitute CA, = A, A; Epo = 0.28 Ey, 


= (0.256) [160 + 70) Btu/sq ft hr 
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Figure 43. Viewing Factor 


= 58.9 Btu/sq ft hr 


0.25 58.9 7 225 
re (Gey - {58:9 *1.25) = 455R = -5F 
0 4 THRU Y eae 2 


(Ref. 13) 
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APPENDIX I 


KEKKEKKEKKKEKEKREKREKKEREKE y r 
* * 
«x WALTERS, W * 
* THESIS RESEARCH * 
* 6 JUNE 1990 ‘ | 
* NONLINEAR SYSTEMS * 
* * 


KRKKKKKEKERKKEEKKREEKEEEK 


** *DURPOSE*** 
THIS PROGRAM COMPUTES DISPLACEMENTS, VELOCITIES,AND 


ACCELERATION 
USING LINEAR ACCELERATION METHOD 


***VARIABLE DEFINITIONS*** 


ta 
} 


A= Acceleration 
= Velocity 
= Displacement 
= Mass 
= Spring constant 
= Damping 
= Force ‘ 
FS= Spring force 
KSTARS= Delta K* 
PSTAR= Delta P* 
DELT= Time interval, delta t 
DELU= Change in displacement 
DELV= Change in velocity 
DELA= Change in acceleration 
DELP= Change in force 
A0= Coefficient of amplitude for force function 
Omega= Frequency for force function 
Al= Coefficient for X**3 term for stiffness 
function 
Bi= Coefficient for X**2 term for stiffness 
function 


Cl= Coefficient for X term for stiffness function : 





c 
Cc 
Cc 
c 
Cc 
c 
Cc 
c 


*DELA, DELP, A0,Omega,Al,B1,C1,D1,A2,B2,C2,FRAC,MASS,P(5000) 


c 


c 


aa 


Qa aaaana 


Dis coefficient for constant term for stiffness 
function 

A2= Coefficient for X**2 term for damping function 

B2= Coefficient for X term for damping function 

C2= Coefficient for constant term for damping 
function 

FRAC= Fraction of linear stiffness term used for 


coeff. of cubic term 


***VARIABLE DECLARATIONS*** 


Real*16 T,A,V,U,M,K,C,FS,KSTAR, PSTAR, DELT, DELU, DELV, 


DIMENSION P(5000) 
Integer I,N 


CHARACTER*7 DUMMY 
CHARACTER NAME*8,FN*8 
***INITIAL DATA*** 


***MAIN PROGRAM*** 


=s= CREATE A FILE FOR CHAOTIC COMPUTAIONS 


OPEN (UNIT=1, FILE='ACC.INP’, STATUS=’OLD’ ) 


WRITE (6,*) ‘INPUT THE FILE NAME’ 
READ (5,’(A)‘) NAME 

FN = NAME 

OPEN (UNIT = 12, FILE = FN, STATUS = 


WRITE(6,*) ‘DO YOU WANT INTERACTIVE INPUT FILE, 


0=YES,1=NO’ 
READ(5,*) IYES 
IF (IYES.GT.0) GO TO 500 


PRINT *, ‘ENTER INITIAL CONDITIONS’ 
PRINT *, ‘INITIAL DISPLACEMENT =’ 
READ *, U 

PRINT *, ‘INITIAL VELOCITY =! 

READ *, V 

Print *, 'F(X)=AG*COS(Omega*T) ’ 
Print *, ‘Enter value for AO’ 

Read *, AO 

Print *, ‘Enter value for Omega’ 
Read *, Omega 


Print *, ‘K(X)=AL*X**3+B1*X**2+C1*X+D1! 


Print *, ‘’Enter value for Al’ 
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509 


550 


Read *, Al 
Print *, 
Read *, Bl 
Print *, 
Read *, Ci 
Print *, 
Read *, Di 


‘Enter value for Bl’ 


"Enter value for Cl’ 


‘Enter value for D1’ 


Print *, 'C(X)=A2*V**2+B2*V+C2" 


‘Enter 


‘Enter 
Read *, C2 
Print *, ‘Enter 
Read *, M 

Print *, ‘Enter 
Read *, DELT 
Print *, ‘Enter 
Read *, N 


GO TO 550 
CONTINUE 
READ(1,;*) 
READ(1,*) 
READ(1,*) 
READ(1,*) 
READ(1,*) 
READ(1,*) 
READ(1,*) 
READ(1,*) 
READ(1,*) M 


WRITE (*,*) AO,Omega,Al,C1,B1,A2,52,C2,M 
"NONLINEAR STIFF TERM Bl IS FRACTION OF LINEAR 


PRINT *, 
TERM D1’ 
PRINT *, 
READ *, FRAC 


Bl= FRAC*D1+D2 


PRINT *, 
PRINT *, 
READ *, U 
PRINT *, 
READ *, V 
Print *, 
Read *, DELT 
Print *, 
Read *, N 
CONTINUE 


value 
value 
value 
value 
value 


value 


for A2‘ 
for B2’ 
for C2’ 
for mags’ 
for deit’ 


for n, number of time steps’ 


‘INPUT FRACTION DESIRED’ 


‘ENTER INITIAL CONDITIONS’ 
‘INITIAL DISPLACEMENT =’ 


‘INITIAL VELOCITY =’ 
‘Enter value for delit’ 


‘Enter value for n, number of time steps’ 
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- —v eas ee we ee ne a ee 
ee NCE ER A LE IE, 


1 om 


pares 


anaaq anaq ama ann anaann aan 


KKK 


6000 
c 
Cc 
20 
c 
c 
c 





DUMME="R8ALTERS’ 

Do 20 i=3, 

Cali Force(£0,Ozeca,=,2,>) 

Cal: Dazp(A2,52,C2,4,C) 

Call SSPRING(A1,31,C2,D1,0,2,X,FS) 
CALL DFORCE(T, DELP, 2,2 

Solve Zor & 
B=(1/4)* (2 (2)-(C*¥)-FS) 

Solve for KSTAR 

KSTAR=K+( (3*C) /DELT)=+( {6*H) /DEL?**2) 
Solve for PSTAR 


PSTAR=DELP+( ( (6*M) /DELT)+(3*C) )*V+( (3*H)+C* (DELT/2) )*A 


Solve for DELV 
DELV=(3/DELT) *DELU- (3*V)-A* (DELT/2) 
Solve for DELA 
ELA=(6/DELT**2) *DELU-( (6/DELT) *V)-(3*A) 
Write output file 


OUTPUT TO THE CHAOTIC FILE 
WRITE(11,6000) DUMMY,T,U,V,A 
FORMAT (A7,F6.3,£11.4,4E12.4) 


Update U,V,T for next iteration 


U=U+DELU 
V=V+DELV 
T=T+DELT 
Continue 


Write(*,*)T,U,V 
CLOSE(UNIT-1) 
CLOSE (UNIT=11) 
Stop 


End 
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KEREREREEEAEEKEKEAEEEEEKEAEEEAEAARAEAEEEEREAAEEKEAKEENRAEEREEAEKEEK 
Subroutine Force(A0,Ozece,T,=,2) 


Rssicns a value to force based on tire . " 
Input arqucent 
30-Coeziicient,Onmeca-fzecuency,T-Tice, i-Counter 
siput arcucent 
P-Force 


AAMNAAAN ARA 


Real+16 20,O=esa,T,? 
iInteaer = 
DIMENSION (5000) 


Define results 


P(I)=A0*COS(Omega*?) 
Write(*,*)1I,2(2) 
Return 

End 


Qa aaan 


4 
' 


KKK KK IKK KKK KEKE IKKE KE KEKE KKK KE KKK KKK KEKKERKKEREKKEKKKK 
Subroutine Damp(A2,32,C2,V,C) \ 


Assigns a value for damping 

Input argument 

A2,32,C2 - Coefficients of damping function 

U-Displacement 

Output argument 

C-Damping { 


Real*16 A2,82,C2,V,C 
Define results 


aq ananaanaaaq ar#an 


C=A2*V**24+B2*V+C2 

Return 

End 
c ‘ 
Cc 


KKREKK KKK KKK EKER KKK KKK KR KK KEK KEKE KK KKK KEE KK KREKEEKRERKKKD KK 


Subroutine FSPRING(A1,B1,C1,D1,U,1,K,FS) 


Calculates the value for FS 


Input arguments 















A — 





U-Disolacexents 

Bi,31,C1,D1- CoeZzicients ior stifiness function 
Z-counter 

Output arcuments 


FS-Force of spring 


aanagaaanaa 


Reai*16 A1,31,C1,D1,U,FS,xX 


integer = 


Cc 
c Define results 
c 
K=A1*U**3431*U**2+C1*U+D1 
FS=K*U 
Retura 
End 
c 
c 
KKKKKKKEKKKKKEKHKKKEEKKEERERELREKKKKKEKEREKEKKEKKEKEKEKKKKKKKKK 
Cc 
Subroutine DFORCE(T,DELP,I,P) 
Cc 
c Defines the value for the change in force 
Cc 
c Input argument 
c T-Time,P-Force 
c Output argument 
c DELP-Change in force 


Real*16 T,DELP,F,P(5000) 
INTEGER I 
DIMENSION P(5000) 


ann 


IF (I.EQ.1) Then 
DELP =0 

ELSE 
DELP=P(1)-P(I-1) 

Endif 

Return 

End 

c 
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Figure Jl. Fractal Dimension. 


a ee Re me me 


C=.177 


NOISNSWId GoVdS ASVHd 


1.53, 


W 
I 
ad 
T 
Vv 
L 
0) 
Vv 
a 
i 


Figure J2,. Fractal Dimension. F=83.35, 
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Figure J3. Fractal Dimension. F=166.7, 9 
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Figure J4. Fractal Dimension. F=208.4, 0=1.53, C=.177 
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Figure J6. Fractal Dimension. F 
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Figure J7. Fractal Dimension. F=83.35, Q 
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Figure 38. Fractal Dimension. F 
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Figure J9. Fractal Dimension. 
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Figure J10. Fractal Dimension. F=41.68, %=1.53, C=.0454 
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Figure Jll. Fractal Dimension. F 
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Figure J12. Fractal Dimension. F=166.8, Q 
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Figure J14. Fractal Dimension. F 
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Figure J16. Fractal Dimension. F=83.35, 9 
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Figure J17. Fractal Dimension. F=166.7, 9 
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Figure J18. Fractal Dimension. F=208.4, 0=13.8, C=.0454 
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